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Extreme correlations arise as the limit of strong correlations, when the local interaction constant 
U goes to infinity. This singular limit transforms canonical Fermions to non canonical Hubbard type 
operators, with a specific graded Lie algebra replacing the standard anticommutators. We are forced 
to deal with a fundamentally different and more complex lattice field theory. We study the t-J 
model, embodying such extreme correlations. We formulate the picture of an extremely correlated 
electron liquid, generalizing the standard Fermi liquid. This quantum liquid breaks no symmetries, 
and has specific signatures in various physical properties, such as the Fermi surface volume and 
the narrowing of electronic bands by spin and density correlation functions. We use Schwinger's 
source field idea to generate equations for the Greens function for the Hubbard operators. A local 
(matrix) scale transformation in the time domain to a quasiparticle Greens function, is found to 
be optimal. This transformation allows us to generate vertex functions that are guaranteed to 
reduce to the bare values for high frequencies, i.e. are "asymptotically free". The quasiparticles 
are fractionally charged objects, and we find an exact Schwinger Dyson equation for their Greens 
function, i.e. the self energy is given explicitly in terms of the singlet and triplet particle hole vertex 
functions. We find a hierarchy of equations for the vertex functions, and further we obtain Ward 
identities so that systematic approximations are feasible. An expansion in terms of the density 
of holes measured from the Mott Hubbard insulating state follows from the nature of the theory. 
A systematic presentation of the formalism is followed by some preliminary explicit calculations. 
We find a d-wave superconducting instability at low T that formally resembles that found in the 
resonating valence bond (RVB) theory, but with a much reduced Tc. 
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I. INTRODUCTION 

We present a theory of an extremely correlated quantum liquid (ECQL) in this report. We believe it to be both 
necessary and useful to make a distinction between the class of systems studied here, and the more commonly addressed 
strongly correlated electron systems (SCES). The latter deals with models such as the Hubbard model for transition 
metals or copper oxide systems, and Kondo or the periodic Anderson models as germane to heavy Fermi systems, 
and also refers occasionally to the t-J class of systems. We use the term ECQL exclusively for systems such as 
the t-J mode]—, having much stronger- even extreme correlations built into them. The origin of the difference is the 
infinite Coulomb repulsion at each site, where double occupancy is prohibited, not just discouraged as in the Hubbard 
model. The result is that we must deal with a genuinely different problem, both physically and mathematically; here 
the fundamental operators are no longer canonical Fermions but rather Hubbard operators. As we describe below, a 
systematic theory of the ECQL can be built, in parallel to the Fermi liquid (FL) theory for the more conventional 
interacting electron systems. The ECQL is described below as a coupled spin and charge liquid that breaks no 
symmetry, while it accommodates the extreme correlations of the t-J class of models. Its Greens functions and 
vertex functions can be defined in a specific way that is different, and fundamentally more complex than in weak 
coupling (FL) models. In this formalism, the instabilities of the ECQL in charge, spin and superconducting channels 
can be studied systematically. The extremely correlated quantum liquid described here, is a strong coupling entity 
and in contrast to the Fermi liquid, it cannot be related in an adiabatic fashion to a free fermi gas^. Although the 
t-J model is obtainable from a non interacting Fermi gas by turning on the J parameter, one is also obliged to turn 
on the Hubbard U term all the way to J7 = oo. Adiabaticity is lost in the passage to the infinite U limit; Appendix D 
illustrates this breakdown of adiabaticity within the context of the atomic limit of the Hubbard model, and outlines 
the general argument for the invalidity of the Luttinger Ward theorem in the extreme correlation limit. The ECQL 
is expected to describe the Physics for a sufficiently large U, and not just U — oo. 

The theory of strongly correlated electronic systems, as opposed to the extremely correlated matter studied here, 
has received considerable attention in recent times. A notable success has been the idea of large dimensionality 
where the dynamical mean field theory^i^ leads to a controlled set of calculations that are useful as well as predictive. 
The t-J model Eq([T]) has also been studied^, but perhaps not so widely as the Hubbard model. The origin of the 
t-J model is the subject of several studies summarized in ReS. We emphasize that it is not merely a descendent 
of the the Hubbard model (upon using a large U expansion), but rather has an independent origin via the down 
folding of multiband systems. Thus t and J may be viewed as independent parameters, rather than being fixed 
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through J = At^/U as in superexchange theory. The possibihty of an expansion of the physical quantities in terms of 
Mott Hubbard holes, i.e. the departure from half filling, has been anticipated in 1986^, but a systematic procedure 
remained undiscovered. Early studies^ii, have struggled with the technical difficulties of the non canonical nature of 
the Fermions. The proposal of Anderson^ in 1987, that the High Tc systems are described by this model, has led to a 
revived and wider interest in this model. Alternate techniques such as the auxiliary field method have been employed 
to deal with the constraints^, but give divergent views of the outcome of extreme correlations. The present study is 
motivated by this conundrum, as well as by the possibility of its applicability to certain systems found in nature. The 
predicted behaviour of the ECQL has certain specific signatures that distinguish it from the FL, so that it should be 
relatively straightforward to verify its applicability to a given system. 

The plan of the paper is as follows. In SectionlUwe present the calculation of the Greens function of the t-J model, 
along with the background definitions of the Hubbard operators and the Schwinger framework. In Section IIIIl the 
exact Schwinger Dyson equation is given after the sources are turned off. In Section HVl we solve for the local Green's 
function. In Section |V] we outline the computation of the vertices, with detailed results given in Appendix [G] In 
Section |Vl] the Ward identities are developed and current vertices defined. In Section IVIll we initiate the programme 
of successive approximations to the Greens function and list several consistent schemes. In Section rVIIII we present 
a calculation of the d-wave superconducting instability of the ECQL, using one of the consistent schemes for the 
Greens functions. In Section HXl we present a detailed picture of the quasiparticles which are argued to be fractionally 
charged. The argument for the changed Fermi surface volume is given further in Appendix |D] We summarize the 
results in Section 1x1 

Appendix |X] and its subsections lA 11 IA2I and Appendix [B] provide the details of the functional derivatives and the 
functions that arise in that process. Appendix Ogives the details of the Nozieres relations that play an important role 
in giving us rotation invariance of the theory. Appendix |E] provides the details of different zero source limits that are 
available. Appendix |F] contains details of the connection between the quasiparticle susceptibilities and the physical 
particle susceptibilities as well as some important sum rules that set the scales for the susceptibilities. Appendix 
Ogives all the vertices computed after throwing out higher order vertices, i.e. Finally Appendix IHl gives the 

conventions used for various Fourier transforms. 



II. THE BARE AND QUASIPARTICLE GREENS FUNCTIONS OF THE t-J MODEL 

We study the t-J model given by 

if = - ^ ^ijC^ ,o-Cn,o- + ^ X! "^idi^r, ■ Sp. - ^nf.np^ + npj - n?,, (1) 

where ct ^ creates a electron at the site i which is Gutzwiller projected^'^ to the subspace of single occupancy, and 
the other symbols have their usual meaning. The Gutzwiller projection creates the major technical challenge in 
this problem. Projected electrons no longer satisfy canonical anti-commutation relations and the power of Wick's 
theorem^i is lost, therefore a Dyson equationJ^ with straightforward expansion in terms of a free Greens function 
is not possible. Lost too are the appealing Feynman diagrams that encapsulate standard FL theory. The projected 
electrons satisfy instead, a set of graded Lie algebraic commutation relations. The latter are compactly expressed 
in the notation of Hubbard in terms of the X operators that are summarized in Section III Al Section III Al contains 
the details of the definitions of the Greens functions, and of the equations of motion calculation, using the powerful 
technique invented by Schwinger, Martin and their schooU^"— . We find that despite the non canonical nature of the 
Hubbard operators, we are able to obtain equations that are at the same level of complexity as those of canonical 
Fermions, with some unavoidable embellishments^^. We therefore conclude that in the case of projected electrons, 
the Schwinger-Dyson expansionii is yet possiblei^. This work studies the resulting equations and their consequences. 

A. Basic framework using the Hubbard Operators 

Let us define the Hubbard operators acting upon each site as projected Fermi operators 

= \a,}{a,\, X°^^ = |0)(ai|, = |ai)(0|, (2) 

in terms of the three possible states at any site |0), | t)i I i)- The doubly occupied site is forbidden, and these operators 
do not connect forbidden states with the allowed ones. Thus an important statement of completeness at any site is 
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the relation 

x^^ = i-Y,xr- (3) 

(7 

The first member X^^'''^ is Bosonic while the other two are Fermionic with respect to their commutation relations 
at different sites. The fundamental anticommutator for the Hubbard operators at different sites is most conveniently 
expressed as 

= 5ij%,..-CT^CT2Xr']. (4) 

Also note the basic commutator 

[X°-.\x;^-^] = c5.,c5.3ai^r^ (5) 
The Hubbard algebra at a given site is defined by 

XfXf = StcXf (6) 

In brief we may visualize these Hubbard operators in terms of familiar (canonical) electronic operators ct ^ creating 
an electron at site i with spin cto, via the non- linear constructs 

X^-"' = {l-nr,,j4^^^^cr„., (7) 

where aa = —da here and throughout this paper. The factors of (1 — rif^ get rid of doubly occupied sites. Although 
one can work with the c's and project doubly occupied states out, it is optimal to work with the X operators. As 
Hubbard pointed out, the manifold of lower Hubbard band states defined by the reduced Hilbert space excluding 
double occupation, is mapped into itself under the action of X's. 
The Hamiltonian is expressed in terms of the X operators by 

The last (trivial) term is a shift of the chemical potential, and is added to make the equations more compact. The 
third (J dependent) term is rewritten as 



^ = i E -^ii^p^i - xr^xp^^} = J E Jijixr^xp^- - xp^^xp^^y, (9) 

ij ij 

where P permutes the spin indices. 



B. Calculation of the Greens function 

The Greens function is defined as 

g„,„M = -(T(xr^(n)xr°(T2))), (10) 

where T is the time ordering symbol, and for an arbitrary operator Q we define 

Thus, ^ is a 2 X 2 matrix in the spin space. The time dependence is given by 

Q(r) = e^"Qe-^". (12) 
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In further work we need to add a source teri»i^ via the operator A 



I dTVp''^T)Xp'^^T), (13) 
Jo 

with the same t dependence of the operators as in Eq (jl2p . and an arbitrary function of time VJ^'^^ (r) at every site. 
The Greens functions are now no longer time translation invariant, and are defined as 

. rr[e-^^r(e-^Xr(n)Xr"(r2))] 

This version of G satisfies the Kubo Martin Schwinger boundary conditions-'^'^ as may be readily verified 

Gcncr2[riO^ ,r20] = -ScritTa [''i/^, '"2O], (15) 

so that one may perform a Fourier series with odd integer Matsubara frequencies only. We will use the convention 
that the Greens function is written as the difference of frequencies after setting the source ^ = 0, whereas for A ^ 
we will display a function of two separate times. 

More generally for any variable we define a modified expectation 

((Q(Ti,r2,..))) = rr [e-f^"T{e-^)] ' ^^^^ 

with a compact notation that includes the time ordering and the exponential factor automatically. Thus 

= (17) 
From this the variation of the Greens function can be found from functional differentiation as 

^^^Jj^((0(r2))) = ((0(r2))) {{Xp'-Hn))) - {{Xp'^^in)Qir2))) (18) 
and we note the important commutator: 

j 3 

We note the similarity in form between the t and J terms above, one can be transformed into the other by hipping 
the spatial indices i,j on the X operators; this symmetry persists in the following equations as well. From this point 
we will use an Einstein type convention, we sum over all internal repeated indices, while leaving the external indices 
fixed. By summing over an index, a spatial sum over the lattice and integration over imaginary time < Tj < /3 is 
implied. The external indices (both space-time and spin) are recognizable since they appear in the LHS of all the 
equations. 

Let us compute the time derivative of the G. For this we need the derivative 

drj(e'^Xt\n)) = -T(e-^[Xt\n),H]) +Vr"- {nW (e-^[Xr'- {n), xf {r,)]) . (20) 
This follows from the definition of the time ordering and the form of A. Using this we find; 

dr^g^^., [^, f] = -Sin - rf)6,j{{{s.,., - '^.'^Z^"'') )) + (([X^' (n), H] X^^^rj))} 

-Vr"^(ri)g.,., [*,/]. (21) 
We further use the abbreviations, 

S[hj]^SijS{Ti-Tj), t[i,j] ^ tij 6{Ti - Tj), 

J[iJ] = J,j 5{n - Tj), V;-'^" = V",^^'[Tr]. (22) 

We next introduce a useful and convenient notion of "k-conjugation" of any matrix M, denoted by M*^, such that 
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This is the time reversal operation, but confined to the spin space indices and excludes transforming the space-time 
indices. The k-conjugation of a matrix in the spin space produces the transpose of its cofactor matrix, and thus the 
inverse of any matrix is proportional to its conjugate, as in Eq (f24| . More explicitly we see that 

M ■ M'= = det M 1 

(My^ = M^ (23) 

^ detM ' ^ ^ 

where 1 is the identity matrix in the 2x2 dimensional spin space. 

We need an object A[i] that plays the role of a dynamical Gutzwiller factor in this theory. It is defined by 

1 

detA\i] 



^^Ll^] - XTaT7T<..W- (24) 



where the second and third hnes follow from Eq (j23p . We also define a matrix functional derivative operator 



DauaAA = f^lf^2^^;^p^. (25) 



In terms of these, we find the equation of motion (EOM) 

{dr, - ^J)Qa,ai[^J] = -S[i, f]{{5^^ - a^a f X^" ^ )) 

+ lJ[^,J] (( {S..^.,~a,a,x;''^) ^^(r,) X^/^ir^)) 
-vr^in) g.^.A^J]- (26) 
We employ a useful relation with an arbitrary operator Q: 

{{{6,^,, - aaa.X:^^-) Q)) = (A,„.Ja] + D.^.Ja]) ((Q)), (27) 
to rewrite Ea (l26p in component form as 

(a.. - m)^..., /] = -5[^, /] A..., W - Vf'"^- (n) [i, f] 

+ t[^,J] { (AW + D[i]) ■ /] L^,^ - { (A[j] + D[j\) ■ g[z, /] }^^,^ . 

(28) 

This may finally be written compactly in matrix form as 

(^n - ^^)Q[^, f] = -A^, nm - H • Q[^, f] - X[i,j] ■ g[j, f] - Y[t,j] ■ g[j, /], (29) 
where we used the definitions 

X[i,j] = -t[t,j]D\i] + ^J[i,k]D[k]6[i,j] 

Y[t,j] = -t[i,j] AH + A[k]5[t,j]. (30) 

The space-time indices are displayed but the spin indices are hidden in the above matrix structure. 

We next perform a scale (or local gauge) transformation with a space time spin dependent factor described below. 
This scale transformation is a key step in our work, and it is important to appreciate its motivation. If we work 
with the EOM Ea ip^ . the resulting vertex, i.e. schematically the object — turns out to have pathological 
"overhangs". By this we mean that the vertex will contain not only g, but also g~^, i.e. the putative "self energy", 
and hence the resulting Schwinger Dyson equation will be ill formed^. The origin of the difficulty is that the coefficient 
of 6[i, f] in the RHS of Ea (P^ involves A[i] which is time dependent, and essentially made up of g. This in turn is a 
manifestation of the non canonical nature of the projected electrons. The factor A[i] is physically understandable as 
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arising from the spectral weight contained in the lower Hubbard band, which is less than unity. On the other hand, 
canonical Greens functions contain a spectral weight of unity. These observations are reflected in the coefficient of ^ 
in the limit of high frequencies of the Greens functions changing from unity in the canonical case to 1 — for the 
t-J model. 

We resolve this difficulty by a local space time dependent scale transformation. With this, we eliminate this 
inconvenient factor of A[i] by a multiplicative process in the time domain. After the removal, we uncover a new 
Greens function G corresponding to effective canonical electrons or quasiparticles (QP) Eq([7]) that lie underneath. 
One could remove this factor A[i] in many equivalent ways, such as symmetrically or from the right, we choose a left 
sided transformation for maximal ease of computation^^. Let us write 

g[ij]^m-G[i,f], (31) 

so that the EOM Eq (P5)) becomes after some rearrangement: 

- 5[i, /] 1 = {(a,. - + + <^m,]] 1 + A-M*] ■ x[i,j] ■ A[j] + A-i w ■ j] • a[j]} .g[j, /] 

a>, = ^-\i]-{drAm- (32) 

Here G is the underlying canonical Greens function (with spectral weight unity), V is the transformed source field, 
and $ arises from the time derivative. Appendix [B] summarizes its properties, and demonstrates that it is negligible 
on turning off the source terms. 



Calculation of the inverse Greens function 



We next turn to the task of finding the inverse Greens function. Using the methodology and results detailed 

in Appendix lA} (especially Eq (p4l lAip ). and with :Q: symbolizing a right (i.e. normal) ordering of the functional 
derivative contained in the matrix operator Q, we may rewrite Ea p2p as 



- 6[i, /] 1 = {dr, - A* + K: + $0 • G[i, /] - t[i,j] A[j] • G[j, /] + i J[z, k] ^i[i, k] ■ AH • G[z, /] 

+ {ii[i,k] ■ v[k,i] + :/i[i,fc] • «]:) • G[iJ] (33) 

Here we denote 

= A-iW-A[j] , (34) 
v[k,]] - ;P[j].A-i[.7];-A[ft] , (35) 

^^Ti.a^M = 0-1^2 (36) 

The matrix products are indicated by the center dots, and the terms under the overline symbol indicate the extent 
of terms over which the derivative acts. Using the results detailed in Append! jXl Ea (jA3|) and Sectior HVl Eg ((62)) . we 
express various objects in terms of G (rather than Q): 



A[j] = ^(l_G'=[j-,j]), where 



7[j] 

7b-] = l-det[G[j-,j]]. (37) 

With these identifications, we have converted the problem to one only involving G and V at this point, and jettisoned 
all reference to the original Greens function Q and the original source V. 
Our next aim is to find an equation for the inverse^^- of G defined through 

G[z,j]-G-i[j,fc] = G~%]]-G[],k]^t5[i,k] 
G.,,aAhj]G-l,,[],k] = J.„,3<5[*,fc]. (38) 



8 



It is useful to define a "susceptibility" type three point object x 



xl:Z[p.<i;r] ^ (39) 



and a vertex function T 



r--[p,<j;r] = BO that 



The relations of these susceptibilities with the physical ones are detailed in Appendix [F] Let us right multiply Eq(|33 
by G~^[f, m] ( and sum over /) so that 

-G'~^[i,m] = {dr, - ii + Vi + '^i) 5[i,m] 1 - t[i,m] {v[i,m] + A[m]) 
+ - J[i, k] {^[i, k] ■ v[k, i] + /^[i, k] ■ A[i]) 5[i, m] 



-t[i,j] \V\l]■^l[^,3i■ ■G[j.J]-G-\f,m] + \jV,k] \ti[i,k] ■ V[k] ■ ii[k,i]: ■ G[i., f] ■ G'^f M 



(41) 



More explicitly, the derivative excludes operating upon the factor G ^ [/, to] , since the overline excludes that term. 
The detailed calculation is presented in Appendix [A) We use Eq (jAl4IAl5P and Eq ljAlOp to rewrite this expression. 
We split the Greens function into the following form with Gq as the bare quasiparticle Greens function and S as the 
quasiparticle self energy, 

G~'^[hj] = G'^^[hi]- S[i,i], where 
-G^%j] EE (a._;, + ^, + i(^j,_,)(l-^))5[^,j]-i[^,j](l_^). (42) 



The self energy is given by 



2' 

k 



77 

S[i,3] = mhj] - t[i,3] {(A[j] - (1 - -)) + 

1 7i 

-t[i,k] Q[i,k,j] + 8[i,j]-J[i,k] {^i[i,k]-^[i\-{l--)} 
+ ^^[«' k] ^[i, k] ■ i^[k,i] + i J[i, k] ^i[i,k] ■ Q[k,i,j]. 



(43) 



Regarding nomenclature, we call the expression Go as the bare- rather than the unperturbed- Greens function, and 
the rest as self energy, since there is no perturbation theory involved here. The result Eq (l42l43p is the exact Schwinger 
Dyson equation including the source field, and is central to the calculation of the vertices. Calculating the vertices 
proceeds by taking functional derivatives of Eq(P5| , and leads to a standard hierarchy of equations for higher order 
verticesi^ii^. We should also clarify that in contrast to the nomenclature of Fermi liquids, the quasiparticles as defined 
here are not infinitely sharp in energy, they decay through relaxation processes that we describe via the imaginary 
part of their (analytically continued) self energy. Thus the quasiparticles of the ECQL have a coherent part as well 
as a incoherent background part that emerges in complete analogy to the electron Greens function in the FL. 



III. SCHWINGER DYSON EQUATION AND THE LIMIT OF VANISHING SOURCES 

Let us gather various objects here when the sources are turned off i.e. V 0. We recover all the symmetries in 
this limit, translation invariance allows us to perform Fourier transforms in space and time, and rotation invariance 
leads to important simplifications of the vertices and susceptibilities. Due to the rotation invariance of the extremely 
correlated quantum liquid, the three triplet states of the particle hole pair are degenerate. This leads to important 
identities that parallel the relations in the Fermi liquid. We denote these as Nozieres relations and discuss them 
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further in Appendix [Cl where the abbreviated superscripts and subscripts of the vertices are defined. Therefore as 
F -J- 0, 

(44) 

and for any object Q = Xi ^ i or F we have a spin decomposition: 

(Qcri,CT2 > X sj rsj r)(i) I A _ n(2)i_i_\ x x _ n(3) 

= Q(i' + Q('^ (45) 

The object Qg refers to the particle hole singlet channel, corresponding to a charge density variable and Qt to the 
particle hole triplet channel, i.e. a spin density variable. 

Let us note that in the limit of vanishing sources, the physical projected electron Greens function Q and the 
quasiparticle Greens function G are simply related as 

g[z-3] = (l--)G[z-j] and 
g[k,iLu„] = (l-|)G[fc,^L^„]. (46) 

The object v is now a function of the difference a — b, and due to the multiplying factors t[a, 6], J[a, 6] in the 
expression for G, we usually need this at equal times Ta — ti,. We Fourier transform using the convention in Appendix 
El as 

:/o-i,crJa, 6] = 6a^^a^i^[a, ^] = i^[q] exp -i{a - b)q 

p 

1 n 

,y[a,b] = —l-{(!^^l)^W[a,a;b] + ^x^^^[a,a-,b]-x^^^[a,a-,b]}. 

{{Sa ■ Sb) + ■^{{naTib) - n"^)}, with ra=Tb 



1 - f " "' 4 



1 n 1 - - 

K?] = - i,)xs[q] - ^ xAq\- (47) 

2 I 2 1 — n 



Further on turning off the sources, we find from Appendix lA 21 

Qai,a2[r, s,m] = aiaafJ-cr^.crd'^, s]Ga^^crA-'^,k]Tgl]gl[k,m-,r] 
S^,,a2G[s,k] {F(2)[fc,m;r] -F(^)[fc,TO;r]} 
= <5,,.,,G[s,A:]F(P)[fc,m,r]. (48) 

As explained in Appendix [Cl the superscript (p) stands for Cooper pairing channel and is a specific linear superposition 
of the singlet and triplet channels F^p^ — ^(F^ — SFf). Hence we may write the EOM as 

~G-Al,j]^{dr^^^i + <P,)6[t,j] 

+ ^J[i,k] (i.[fc,*] + (l-|)) S[i,j] + ^J[i,k] G[i,l]r^P\l,r,k]). 

(49) 

Using the convention for Fourier transforms in Appendix [Hi and with the Fourier space version of Gq Eq(I 



G^^k] = iLOu + /i - (1 - |)(efe + ^ Jo). (50) 
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we write the exact Schwinger Dyson equation for the Greens function 

G-^[k] - G^^[k]-S[k], 

S[k] = ^ + + + G[q]{T,[q,k]~iVt[q,k]}- 



2 ■ 
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(51) 



The self energy exphcitly contains the electron dispersion Sq or tij as a multiplicative factor here. This is quite unlike 
the standard case of the electron liquid, where the entire dependence on the electronic dispersion is hidden in the 
form of the bare propagator Go , and the self energy is a functional of only the (full) Greens function^-^i^Sr— . 
In order to evaluate the first term of S in Eg (1511) . we need a few definitions at this point. Let us denote 

Xs[Q] ^^Xs[q,q + Q], Xioc ^''[^l = ~f {l_ ' 

- — E ^P+Q ^«[^]' Jp^—J2 Jp+Q Xs[Q]- (52) 

Xloc Q Xloo Q 

Similar definitions hold in the triplet channel. In appendix IF] we show that the local object Xioc is identical for both 
triplet and singlet channels. From Eas (l471 1511 and [5^ . and further using the definition of the renormalized energy 
type variables 



Ek = Sq+k iy[q] = - 



3 71 J 1 71,(1 — n) 



4(1-1) 4 (1-f) 



el 



we find 



^ 1 T ri 3 n . ln(l-n) 

'^k = ^Y. - i (T^l) -^fe + i (1 „ I) ^fe (53) 

5W = (i?fe + ijo) + + Ijk-q) G[q] { Ts[q,k]-3rt[q,k]}. 
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(54) 



By construction, the variables Jp, Jp,ep,ep are weighted averages of the bare objects Jp,£p with momentum depen- 
dent but static weight factors, and hence it is reasonable to view them as correlation adjusted versions of the bare 
dispersions. The precise relation between the bare dispersions and Ek-, Jk is particularly simple when the hopping is 
only nearest neighbour; it involves the the spin spin and density density correlation functions at nearest neighbour 
distances, and is given explicitly below in Ea (j58p . Specializing to only nearest neighbour hoppings, we see from 
symmetry that the form of the band dispersion remains a simple tight binding one. Therefore we write a convenient 
notation of various objects for the simple cubic lattices: 



— E<^°^'?2; Xt, s\q\-, (55) 



The parameters Xs and Xt can be found in terms of equal time correlations by carrying out the frequency and spatial 
sums, and using Eas (|F5IF8|) . we find 

Xs = — - — "-^—^ , and Xt = - — {Sp^ ■ Sf^+fj), (56) 

n — 6 n 

where 7f is a nearest neighbour vector. In terms of these, we may write 

£\ = Xtek, £l=£sJk, Jl=XtJk, Jk=XsJk. (57) 

Thus we can express 

Ek = {(-%•-%) + t(("o ^i^) -^^)}' 

Jk = {( ^0 ■ Sfi) + ^(("0 n^) - n^)}. (58) 
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where ff is a nearest neighbour vector and the correlations are at equal times. Thus the nearest neighbour charge 
and spin correlations directly influence the effective band width. The spin correlation has its largest magnitude 
near half filling, and can be positive or negative depending upon whether ferromagnetic, i.e. Nagaoka Thouless type 
correlations, or the more usual antiferromagnetic correlations prevail. Near half filling the density term in Eg (1581) is 
suppressed, while the spin correlation term survives and gives the dominant contribution. This correction term to 
the energy dispersion E is of the same form as e but has the opposite sign and hence leads to an important band 
narrowing via the correlation functions indicated in Eg (1581) in this theory. 

Later we will find equations for the susceptibilities, so that these correlation functions can be found in terms of the 
vertices and the Greens functions self consistently. Successive approximations for the vertices will be formulated, with 
each approximation providing a complete calculational scheme. We will find that the susceptibilities Xs.t[Q] vanish 
generically as 1 — n near half filling, in a similar fashion as Xio^ and hence these renormalized energies tend to a finite 
non zero limit at half filling. 

Eg ([54)) is the exact Schwinger Dyson equation for the t-J model. This breakup of the inverse G into Go and a 
self energy type object S is to some extent arbitrary, since there is no a priori notion of an unperturbed Greens 
function. Our breakup has the natural advantage that the resulting bare vector vertices are frequency independent 
and in agreement with independent arguments (see next section where the bare current vertex is obtained in Eg (jl03p ). 
Further they satisfy the Ward identities as noted below Eg (l97t . The S object, after the standard analytic continuation 
to the retarded self energy, provides a correction to the quasiparticle energies at the poles of G through its real part 
and more importantly to the decay of the quasiparticles through its imaginary part. The self energy in our formulation 
is developed in terms of the Greens function G and is explicitly a functional of G. However, and in contrast to the 
usual weakly interacting Fermi system formulatioit22^i2ii^, it is not a universal functional, i.e. it does depend upon 
the bare dispersion tij , regardless of how we define the self energy. 

Finally it is important to note that this Schwinger Dyson equation is well formed in the sense that the vertices F 
are guaranteed to be well behaved (i.e. finite) for u oo, since we have avoided the linear oj dependence of F that 
results, if one does not extract the time dependent factor A[i] as in Eg ((3T|) . The equations for the vertices are derived 
and discussed below. 



IV. SOLUTION OF THE LOCAL GREENS FUNCTION 



We next consider the various local Greens functions and their interrelations. Start from 

g[j-,j] = A[j]-G[r,jl with 
A[j] = t-g''[j-,j], so that 

m = 1-G'[r,j]-A'[j] (59) 

This equation is easy to solve when we iterate once, 

A[j] - 1 - G'^b-, j] + G'=[r , j] • A[j] • G[r,jl (60) 

we see that the inhomogeneous term 1 — & commutes with A iteratively, so that we may as well rewrite the second 
term in the RHS as A[j] ■ G''[j~,j] ■ G[j~,j] = det[G]A[j], whereby the solution is noted in Eg(|37)): 

A[j] = A^(l-G^[r,j]), where 
7b] 

7[j] = l-det[G[r,j]]. (61) 



We also require the inverse 



A-Mj] = ;^(l-G[r,j]), where 



7lb] = det[l- G[r,j]]- (62) 

These local objects are the dynamical analogs of the Gutzwiller projection factors^*^, and we note their physical mean- 
ing in terms of the bare (quasiparticle) charges n[i] {nQp[i]) and spin densities S[i] {SQp[i]) defined in Egs (ll44ll5ip 



(1 - inH)2 - S[i] ■ 



(1 - ^n[i])^ ~ S[i] ■ S[i] ^ 
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When we turn off the sources, the ECQL state has the following local Greens functions expressed in terms of the 
number density at each site per spin n — N/{2Nsites)'- 



77 — T) 

e[r,j]-2 G[r,j] = Y^ A[j] = {!--), 



7b] = 7lb] - (63) 



The vanishing near half filling of 7 and 7I, as (1 — n) and (1 — n)^ has important consequences, and leads to the hole 
density expansion reported below in Section IVII El 

V. CALCULATION OF THE VERTICES 

The vertices can be obtained from the general equations for the Greens function Eq (|^ as 

rW[i,j;fc] = 5[i,J]S[i,k]+^-^Sn[^,J]] 



= ^W[^,J•;fc]-^(2)[^,J;A:]. (64) 
Thus the bare vertices are simple, writing them in Fourier space we find 

rs[pi,P2] 1, Yt\v\-,V2\ 1- (65) 



The vertex corrections arise from the self energy part of Eq(|64p. and contain several terms. We provide the results 
of the long calculation in Appendix [Gj There all terms are retained, with the proviso that the higher order vertex, 
obtained by differentiating the three point vertex are set to zero, i.e. — t- 0. 

While this complete solution is available in the Appendix [Gl it is crucial to understand the relative order of 
various terms, in order to make sensible approximations. We next provide a set of calculations that give us such an 
understanding, we evaluate the derivatives of the basic elements A[i], and occurring in the self energy, 

and show that these have an explicit dependence upon the hole density b — 1 — n, measured from half filling. The 
terms can in fact be grouped in a formal expansion in the inverse hole density A = j^;;;. In addition there is an 
implicit dependence on n in all terms, and certain terms vanish with b. Thus the final result for the vertex obtained 
by taking a product of the apparently divergent terms and the coefficients thereof either vanish at half filling n = 1 
or remain finite. This gives us an organizing principle for grouping the contributions, the leading terms near half 
filling consist of terms that remain finite at n = 1, and one can throw out the rest, thereby giving us a hole density 
expansion that has been conjectured earlier^. 

Therefore taking the derivatives of the Schwinger Dyson equation with sources Eq (|43p . we find that the expansion of 
the vertices can be arranged as a formal power series in the inverse hole density A, in the form T — l + ro + Ari + A2r2. 
It is implied that the coefficients V j contain terms that are either vanishing at half filling or are finite, and rules for 
recognizing this implicit dependence are provided later. In Appendix [Gl we list the vertices ro,ri,r2. 

The existence of this expansion is fortunate since the most interesting physical regime for a doped Mott insulator 
is in the limit n ^ 1. The factors of A arise from the dynamical projection factors A[z], j] and j] in the local 
Greens functions and the self energy. We next list the leading behaviour of the functional derivatives of the various 
matrix functions ^Ji[^^j]^ j], A[i] w.r.t. the source terms, so that results for the vertices can be assembled together. 
For this section, we will use spin diagonal sources Vm so that the matrices can all be taken as diagonal in spin space. 
We will need to use the Nozieres relations x*^^^ ± x^^'' = Xs. t, in order to obtain the complete the set of derivatives 
in the singlet and triplet channels. After providing the results, we will discuss their relative magnitudes, and also the 
implicit density dependence of the coefficients. 

Let us begin with the less singular terms A and /i, where the functional derivative is explicitly 0(|), and then 
progress to v which is O(^). At some places, we will use the symbol ^ to indicate that the density n has been set 
at unity in all terms except the singular ^.^J^^^ type terms. 
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1. A[i] and its derivatives 



Let us recall for spin diagonal sources f Appendix IE)) 



v^Q 2 — n 



(66) 



Therefore 



nW 



SA^^ [i] 



sv, 



v^o (1 — n) 
y^o (1 — n) 



z ,i:m 



(67) 



^. fc] and its derivatives 

Lets us recall fc] = • A[k], and it may be expressed as 

^in[hk]=g2[k,^] [l - G^^[i~ - G^^[k' ,k]), 
where (?2[«, j] and a related object gi[i,j] needed in the next section are given by 



7b] 



7l[j] 7[i]2 

2 ' 



(1-1)4 



y-i-o (1 — n)3 



7l[j] 7[i] 

n 
2 



(1-1)2 



i^^o (1 — ny 



where we used the definitions (recall Eas (l62l |37)) ). 
We first note the derivatives of g 



^ ^ m 

^g2[hj] 

^ ' 7-n 



v^o (1 - n)4 
v^o (1 — n)3 



+ {1 - -)xs[jj;m] 
-X4z ,i;toJ + (1 - -)Xs[j, j;H 



1 



64(54 
1 

16^ 



[2xs[z ,i;m]+Xs[j,j;m]] 
[Xs[i^,i;m] +Xs[jJ;m]] ■ 



(68) 

(69) 
(70) 



(71) 
(72) 



Notice that the derivatives of gi and g2 are independent of the spin index a. We next compute the derivatives of 
fi[i,k] with spin diagonal sources, by using the above equations and the derivatives of the second and third factors of 
Eq| 



H^^[i,k] 



n 



SVr 



(i-t) 

v^o 2(1 — n) 
1 



[Xt[k ,k;m] ~ xt[i ,i;m] + {l-n){xs[i ,i;in]-Xs[k ,k;m]}] 



— ( Xt [fc ,k;m]~ xt[i , «; m] ) , and similarly 



SVr 



(1-f) 



[Xt[i ,i;m]-xt[k , k;m] + {I - n){xs[i ,i;m]-Xs[k ,k;m]}] 



y-i-o 2(1 — n) 

~ ^ ( Xt[i^,i;m] ~xt[k^,k;m] ) 



(73) 



(74) 



3. v[i,k] and its derivatives 



Using the functions gi 92 in Eq (j69| . We first write i' as 



i^[i,k] = iyi[i,k] + i^2[i,k] 

Uilijk] = fc] fc], and iy2[hk] — g2[i,k] iy[i,k], 



(75) 
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and assuming spin diagonal sources, we find from Eq (|A10|) : 

= -{l-Gn[k-,k]) x<i)r,z;fc]-(l-G;;[fc-,fc]) x^'^i' k]. (76) 

In order to obtain the leading order terms in inverse hole density, we note that the factors 1 — G generate factors 
of 1 — n which lower the order of the term, unless we differentiate these terms first. Hence to leading order, the 
calculation is simply done by ignoring the source dependence of the x's in Eq ()76p . and differentiating the factors of 
1 — G. We thus obtain the leading derivatives of the two objects i^i, 1^2 as follows 



SVl^ Iv^o 1652 

Sui^^ [i, k] \ _ 1 

(5i^2tt [h k] I 1 

SVl^ Iv-^o " 1652 

Su2ff [i, k] I 1 



Xsli" ,r,k] {xs[i',i;m] -Xt[i~,i;m] + xt[k' , k; m]) 

{ Xsli^,i;k] {xt[k~ , k; m] - Xs[i' , i; m]) -xt[i~,i;k] {3xs[i~ ,i;m] + xt[k~ , k;m]) } 
{ Xt[i~,i;k] {xt[k' ,k;m] - 3xs [i~ , i; m]) -Xs[i~,i;k] {xs[i~,i;m] + xt[k~ , k;m]) } (77) 



SVrk^ Iv^o 1652 

Using the Nozieres relations we can combine these and write the final answers in the singlet and triplet channels 



6^ [i, k] 
Sf [i, k] 



singlet 3 

= --^Xs[i^ ,i;m]xt[i~ ,i]k] 



triplet 1 

^ ^{ Xs[i^,i;k]xt[i^,i;m]~xt[i^,i;k]xt[k~,k-,m]} . (78) 



Let us now consider expressions Eq ((67l [74177^ . where we find terms with A, A'^ explicitly appearing ( where A = |). 
The final answer of each term is either finite or vanishes, and to see this we must recognize the implicit dependence on 
6 of the coefficients. The two point susceptibilities Xs, t[hi~',j], with arbitrary arguments are seen to vanish linearly 
as (5 — i> since these involve a particle and a hole in the quasiparticle band, which becomes completely filled at half 
filling. By the same token all three point susceptibilities also vanish linearly as 5 — >■ 0. We should also remember 
that the vertices themselves are non vanishing as (5 — > 0. This behaviour is readily confirmed by a few low order 
calculations. As a useful example of this organization, we next write all the terms in the vertices that have an explicit 
quadratic dependence on A. These are obtained by assembling the explicitly O(A^) terms in Eg (p7[ 1741771) . giving us 

3 1 

^s[i,j;m] = S[i,m]S[j,m] + g {t[i,j] Xs[i,i;m] Xt[i,i;j] - 5[i,j] -J[i,k]xs[k,k;m\ Xt[k,k;i]} 

rt[i,j\m] = S[i,m]S[j,m]~^ t[i,j] {xs[i,i; j] Xt[hi;m] ~ Xt[hi;j]xt[j,jim]} 

"T^ '^t*''^] {Xt[k,k;m]{3 Xt[k,k;i] - Xs[k,k;i]) - 2 xt[i,i;m] Xt[k,k;i]} 

16 

(79) 

The vertex corrections are clearly non zero as 5 ^ since the factor A^ is compensated by two vanishing susceptibilities. 
These can be recognized from the detailed list in the Appendi;fVlEqs (|G2|G3P as arising from r[i,j;m]2 and r[i, j; mjs 
in both singlet and triplet channels. We will see next that these vertices satisfy current conservation and hence are 
studied further in Se dVII El We emphasize that these do not exhaust the set of terms that are finite as (5 0, and 
constitute a convenient subset of the surviving terms. 

VI. ELECTRICAL CONDUCTIVITY, CURRENT CONSERVATION AND WARD IDENTITIES 

We next set up the calculation of the electrical conductivity for the ECQL. We are interested in deriving a useful 
Kubo type expression for conductivity in the t-J model, in terms of the appropriate Greens functions developed here. 
We will also establish Ward identities^"— that relate the current and charge vertices, and find explicit expressions for 
the current vertex. These Ward identities are important relations since they constrain the possible approximations one 
makes, and are necessary satisfy gauge and rotation invariance of the final results. The charge and number densities 
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and currents are trivially related by a factor of Qe the electronic charge, and so we will work with the number densities 
only. The number density ?i[r r] — XZ''[t) has an associated current density given by 

•^>^)^^E'^>+^'?)' J{r+lri)^z^rf{X--^^^X",--XfX$l^) (80) 

v 

(81) 

where ff is the set of unit vectors connecting a site to its neighbours, is implied in all terms. Similar expressions 
can be written down also for the spin density conservation, but since the t-J model involves terms that flip the spin 
explicitly, the form of the resulting equations differ from the charge Ward identities^ . We will be satisfied with using 
the Nozieres relations to verify compliance with rotation invariance; these are necessary and most often sufficient 
conditions for giving us an invariant theory. The lattice version of divergence of these currents is written in terms of 
the densities as 

(v.^ • j(f)) EE z TM - E ^ • '^(''+ ^'?)' = E (K^'x^fU - K^U^'n- m 

The total current operator can be written 

J^Y. J^"-^ =^Y^^Vtf, X^l^X^.- , (83) 

r r.ff 

from Eq ([80|) . and in case of external magnetic fields, we add a suitable Peierls phase factor to the currents. This 
expression for J is necessary for the calculation of the frequency dependent electrical conductivity tenso r^^i^^ of the 
t-J model with r2„ = 2TrkBTn 



n^piinn) = ql / dr e^""^ ( T,J"(r) ), (84) 
Jo 



where we set the lattice constant ao and fi to unity. Here Tq,^ is the diamagnetic part of the response related to 
the plasma sum ruleSS, and may be evaluated as an equal time correlation. As explained in standard texts^ii^, this 
expression should be analytically continued as i^n uj + iO^, in order to obtain the retarded physical conductivity. 
Our task is now to express Ti{iQ) in terms of the Greens functions in a manner that satisfies current conservation. 

Since we are primarily interested in electromagnetic response, e.g. the optical conductivity, the condition Qa^ <^ 1 
is satisfied. Hence we need long wave length response where the lattice structure is not crucial for the current 
conservation laws. However, the calculation is most effectively performed by first noticing that one can find a lattice 
version of the current conservation law Eq (|55)) . valid at all length scales as detailed below. This more powerful 
conservation law^^ reduces to standard electrical conservation law dtp{r. t) + V • J(r, i) at long wave lengths. Our 
strategy is to derive the exact consequences of the lattice version of the conservation law, and then to take the long 
wave length limit to obtain the electromagnetic vertices. 

Using the Heisenberg equations of motion in imaginary time, it is easy to establish the exact conservation laws of 
number density ^^-^[r r] = «(Vp- J{r)) compactly in operator form: 

d^n[rT] = r[f r], (85) 

where T is defined above in Eg ([5^ . The gauge invariance of the t-J Hamiltonian is ultimately responsible for these 
conservation laws. In order to implement the conservation laws Ea (j85p . we first establish the equation for the Greens 
functions as 

n[m])) - ((xo^'X^" rH» = {5[^M - 6[fM)Q.,.AhI]. (86) 

where the terms on the RHS arise from the discontinuities of time ordering implied in the Greens functions^^. The 
two Greens functions in Eq (l86p can be obtained from the action in Eq (|13p by adding terms 



A ^ A + / dT„i {u[m] n[m] + v[m] T[m]) . 



(87) 
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where it can be seen that the effect of u[m] is precisely that of after summing over a and thus gives the 
singlet susceptibilities and vertices. Our next task is to rewrite Ea (l86l) in terms of G. Towards this we recall that 
Q[i, f] = A[i] ■ G[i, f], and note that 



du[m\ dv[m\ J 

on using the conservation law Ea (j85p . Putting these together, we may thus rewrite Eq(|86p in matrix form as 



Su[m\ Sv[m] 

(88) 

where the disconnected parts cancel identically upon using the conservation law Ea (j85p . The second form is specially 
useful for checking the compliance of conservation laws, it is the most compact statement of the Ward identity. 
We define new "7"- susceptibilities" and vertices 

S 5 
x[[i,j;m]^--^G[i,j], rl[i,j;m] = -——rG~^[i,j]. (89) 

The subscripts s are to denote the singlet, i.e. charge nature of these currents. Below we obtain the electromagnetic 
vertices from these by a limiting process. These objects satisfy the real space Ward identities 

{dTr„Xs[i,f;m]-x[['iJ;m]} = {S[i,m] - S[f,m]) G[i, f], or 
{dr^rs[iJ;m]-Tj[zJ;m]} = iS[f,m]-S[i,m])G-%f]. 

(90) 

Fourier transforming these we obtain the momentum space Ward identities 

i{ujp, - WpJ Xs[pi,P2] - xT[Pi,P2] = G[P2] - G[pi] 

iiup,-LOp,)Ts[pi,P2]-rJ[pi,p2] - G-Mpi]-G-1[P2]. (91) 

We next evaluate the T vertex Eq (l89|) . The effect of v[m] is most easily seen as a modification of the kinetic energy, 
since the T operator in Ea (|82p resembles the kinetic energy closely. In fact we can write 

J2v[m]T[m] ^Y.^„Mm] - v[k])X^'X','^ , (92) 

m m,k 

and recalling the definition of the Hamiltonian Eq(IT]) , we can see that the effect of adding this term in the relevant 
time domain equations is to replace t[m, k] — >■ t[m, k](l + v[k] —v[m\). With this observation, we may write the Greens 
function equation Ea (|42l43p with these added sources. We again split the Greens function as in Ea (|^ 

G~'^[i,j] = Go^[i,j]- S[i,j], where 
-Go%j] ^ idr^-f^ + V, + u[i] + ^J[z,k]il-^Mi,j] 

-t[^,J]il + v[J]-v[m-^), (93) 

and the self energy is given by 

71 

S[i,j] = ^,6[i,j]-t[i,j]il + v[j]-v[2]) {A[j]+,y[i,j] - (1 - -)) 
-t[t, + v[l] - v[i]) Q[i, I, j] + ^J[i, k] k] ■ Q[k, i, j] 



1 / n \ 

+ -J[i,k] i^fi[i,k]-,^[k,i]+ ^i[i,k]- A[i] -{!--) j S[i,j]. 



(94) 
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The charge and current vertices follow from taking the derivatives of the above equation. All vertices are naturally 
split into bare and correction terms 

r = 7 + f, (95) 
where 7 = -^G(7^ and f = -^S. It follows from EqdM]) that 

ls[i,j;m] = d[i,m] 6[j,m], 



■jj[i,j;ni\ = -t[iJ]{S[j,m]-S[i,Tn]){l-'^). 



(96) 



In Fourier space we find 



ls[Pl,P2\ = 1, 

2' 



lJ[pi,P2] = (epi -ep.)(l-^). (97) 



We observe that the bare vertices with Ea ((50|) satisfy the Ward identity Eq (|9T|) . 

In order to make contact with our calculation of the Schwinger Dyson vertices, we consider the susceptibility 

T^[i,jM ^ -l-^g„„[^^J] = {{X^-xfT[m]))-{{X^^Xf)) ((TH)) 



g-»Pi(»-m)-»P2(m-j) ^_. (1 _ e'^<^''P^^){{{X°''XfX''^^X°^^^)) - disc}, 



(98) 



where "disc" stands for the disconnected part and while a is fixed, a' is summed over. In the second line onwards 
we turn off the sources, hence we may ignore the "disc" pieces due to parity. We should also note that 77 is a vector 
without a time component, so that the two factors Xj^X^^_^_^ are at the same time t^. We next use the definition of 
the current from Ea ([55|) . and take a specific limit: 

-^J[P ^P, P' + fi„]/p'^p = 



dp'"' 

E f f^din-^rn) (i(T„-T,)e^"-(^'-^^)+'^"(^'"-^^)-*^^-(^"'-^~^)((XO-(r,)Xf (T,)J"(r™))). 



(99) 



Convoluting with the velocity 



^P=^^p-'T.^^*^^^~''''^ (100) 



we see that 



U^^iin„), (101) 



where the factor of 2 in front is from the spin summation. We note that the response functions T can be related to 
the derivatives of G in a straightforward way, as in Appendix |F] We obtain 

'^T[Pl,P2]+'^T[P2-Pl]G[p2] = {^-^) X[[P1,P2] 

Tlb^-pi] = il-!ly^J[p,-p,]. (102) 
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For notational convenience, in the following we write jl = p + Q, with Q = {Q, iJ7„}. Thus we can rewrite Ea (|101l) as 
U^p{^^^„) = 2(1-|) ^e-''^^^ (^_l.|^r[p^p^g]_(l_|)^r[Q]G[p + Q]|^^^ 

(103) 

where the bare current vertex (1 — ■|)fp is a consequence of the definition of the current Eg ([831 . We make the implicit 
assumption that Q in the vector vertex T"[p,p + Q] always satisfies Qqq ^ 1, and all term of 0{{Q°')^) are thrown 
out, i.e. we are in the electromagnetic regime. Observing from the bare vertex Eq(|97pthat jj[p,p + Q], and hence 
the full vertex Tj[p,p + Q] vanishes as Q — 0, and its derivative is odd under parity, we conclude that the second 
term in Eq p03p can be dropped. We may then define the vector (electromagnetic) vertex (with Qao <C 1) 

C[p,p + Q] = -^Tj[p,p + Q] 

X'^[P,P + Q] = G[p]r'^[p,p + Q]G[p + Q], (104) 

with the bare vector vertex 7 as 



d 



is[p,P + Q] ^ ~[i^iT[p,P + Q]] =^v(i-^)> (105) 

so that 



na^(Q) = - 2 ^ e^"-"^ 7f[p,P + Q] G{P\ r>,P + Q] G\P + Q\- (106) 

V 

The Ward identity Ea (j91l) relevant for electromagnetic response now reads 

i^nVs [p, p + Q] - Q • [p, p^Q\ = G-\p^Q\-G-\p\. (107) 

In practice, this infinitesimal form of the Ward identity Eq (|107p is more easy to implement and also sufficient for the 
gauge invariance of the final theory, as compared with the finite form Eq ([9T|l . 

Returning to the Eg (1^5]) , the equations for the vertex corrections F follow from taking derivatives of the self energy 
Ea ((94l) . We may write the charge vertex simply as 

^^\'^rM = [t^^SoA^.j]] = E (^S„^\^,3\\ , (108) 



with the spin index a fixed at say '\. Here the current type sources v[m] are irrelevant and could be set equal to zero 
right away in Ea ([M)) . The second identity in Eq (|108p is a consequence of the rotation invariance relations of Nozieres 
as discussed in Appendix 

In order to obtain the current type vertex corrections, we get two sets of contributions from Eq (|94p . We write with 
say a fixed a, let us say a =t; 

rj,[t,j;m] = {S[t,m]~d[j,m])t[i,j]iy,4i,j] + {S[l,m]-S[i,m])t[zJ]G[l,k] T^P\k,j;z]. 

(109) 

The set of terms r^(,[i, j; to] arise from the explicit factors of v[m] in Eq (|94p . Here S' represents all the terms of 
Ea (p4|) . but with the explicit factors of v[m] omitted. Thus the term Tj ^[i, j;m] arises from differentiating the rest 
of the v[m\ dependence, and originate in exactly the same terms that contribute to Ea (|108p . The only difference is 
that after taking the derivative, a current vertex replaces the charge vertex. Therefore one can obtain the equations 
satisfied by Tj' to] from those obeyed by Ts[i,j', to], in a simple way. In the RHS, we simply replace terms indexed 
by the external symbol to, e.g. rs[p,q;m] by rj[p, q;m] for any p,q. The same separation is also carried out for 
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the electromagnetic vertex F", and may be obtained by taking the hmit as in Eq (ll04p . We note the vertices Tg b in 
Fourier space (with a fixed a) 

q q 

fs b[p,P + Q] = Y.'^p+i ^-M + E ^9 r(^')[g,p + Q]. (110) 
q q 

We summarize the requirements for a "conserving approximation"^'' for the ECQL, it consists of a suitable approx- 
imation for the self energy Eq (|51l54p . the density vertex Ea (|108p and the current vertex Ea (|109lll0p . These must be 
consistent, i.e. satisfy the Ward identity Eq (|107p so that the response is gauge invariant. 



VII. GENERAL PROPERTIES OF THE GREENS FUNCTION AND SUCCESSIVE 

APPROXIMATIONS 



We next turn to computing the Greens functions from the above theory, after making a series of approximations. 
Let us summarize the various objects of interest. Our interest is initially in the positive definite spectral function 
Pg[Qi s-nd the occupation of the k state ruk defined through the representation 

G[k,iuJk]^ ! dv P^^^'''\ pG[k,v] = --'^mG[k,y + i{)+] (111) 

where the positive definite spectral density has a representation 

PG{kA = T^ EK"!^''"^!/^)!' ^Po.+Pp)^{y^ep~e^), (112) 

with Pet — Gxp(— /?£„), Ea is an eigenvalue of the grand canonical Hamiltonian H — pN and Zqc the grand 
canonical partition function. We will define the quasiparticle occupation number toj. through 

mfc = ^( X-0(fc)X"-(fc) ) - ksTY, e-^o^G[A:,zwfc] = j dv f{v) pdkM. f{v) = (113) 



Here a is either up or down, and < nik < 1 fixes the total number of electrons through the sum rule 

1 \ ^ n 

k 



We expect that at T = 0, a sharp Fermi surface of the ECQL exists generically, and is determined by one of several 
criteria in analogy with that of the Fermi liquid. Let us list them separately now (i) The particle number sum rule 
Ea l|114p is one of them, (ii) A jump in toj, is expected at the Fermi surface so that the locus of jumps in this defines 
the "Migdal" version of the Fermi surface, (iii) Another is the analog of the Luttinger Ward sum rule that says that 
the Fermi surface is the domain satisfying the condition ReG{k,0) > 0, or from the spectral representation Eq (|112p 
we require 

5]|(a|A"-(fc)|/3)|^ >0. (115) 

(iv) The final criterion requires that the quasiparticle life time is infinite at the FS, so that the locus of points where 
PE(fc, 0) = defines the Fermi surface. The approximate solutions for G can be tested with these criteria, and we 
shall discuss only the particle number sum rule Ea (jll4l) initially. 

Below we also need the "bubble susceptibility" object xo[Q] defined through 



Xo 



[g] -^GMG[g + 0], (116) 
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which can be evaluated in terms of the spectral functions readily as 

Xo[Q] = -^V/ / rfi^i PG[q, Vl\PG[q + Q, ^2] ^^'''^ ~ ^^^^ ■ (117) 

q 

We next describe a few systematic and consistent (i.e. Ward identity compliant) approximations that can be 
implemented. Detailed numerical calculations within these schemes are currently under way, and the results will be 
presented later—. 



A. Atomic Limit 



The atomic limit is defined by switching off t and J, and so it rather trivial. We make sure that we recover 
the exact answer known in this limit^i^. We set rs[pi,p2] = rt[pi,p2] = 1 and hence pG[k,v\ — 5{v + /x), with 
G^^[fc,iajfc] = iujk + /i, so that the chemical potential is given by 

f{-p) 



2-n 

We compute the susceptibilities from Xs — Xt = Xo with 

71 ( 1 — 72) 

Xo[Q] = -/3<5q,0 feq.O „ II ■)2 ■ 

which agrees with the sum rule in Ea ljFlip . 



B. First approximation 

The first approximation consists of choosing Go [A:] Eq ((50l) and the bare vertices Ea ((97l) . These clearly satisfy the 
full set of Ward identities, In fact this approximation, with J ^ coincides with Hubbard's approximate solution, the 
so called Hubbard-I solution of the equations of motion^S in the limit as U ^ (yo. The particle number can be fixed 
using Ea (|114l) . and we find that the Fermi surface volume encloses a fraction of the first Brillouin zone, rather 
than the Luttinger Ward fraction of j. The quasiparticle number is greater than the bare particle number and this 
feature persists in all subsequent approximations. This violates the Luttinger Ward^'^ volume theoremSS.^ we discuss 
its implications in greater depth later in Section HxF ^ . 

The susceptibility is just the Lindhard function xo — ~XLind[Q], 

1 ^ f{E!°'> -fi)- f{Ei% - P) 

with Eq^'' = (1 — ^)eq. This function vanishes as we approach half filling n — ?► 1 and has the van Hove fingerprints of 
the above defined "large Fermi surface" . It is straightforward to show that this also satisfies the sum rule Ea (|F10p , 

and from Eqs (jF6IF5[) . the physical spin susceptibility is 2/ig Xl[Q] and the physical charge susceptibility or 

compressibility is {1~^)'^xl[Q]- Owing to the structure of the prefactors, the spin susceptibility interpolates smoothly 
between the Pauli and Curie susceptibilities on pushing the density towards half filling n ^ 1, while the charge 
susceptibility vanishes near half filling. 



C. Hartree Approximation (HA) 

The next approximation we make consists of choosing a frequency independent vertex function, and the self energy 
that is also frequency independent. We write 

GH'[k]^iiu;k + fi)-El^\ (119) 
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Specializing to only nearest neighbour hoppings, we see from symmetry that the form of the band dispersion remains 
a simple tight binding one. Therefore for simple cubic lattices using Ea (|55l56l57l58p we write 

^« = {(1 - + ( ^5 ■ S^) + J(K- - or ^0 -(1 - f ) + J^l^^t + (1 - (120) 

The vertices are taken to be the lowest consistent ones 

^T[PuP2] = (1 - ^)(£pi - ep2) + J2^^i+Pi - Sq+P^Hl]- (121) 

g 

We can check that these variables satisfy the Ward identities Ea ipTj) exactly. The object E^f^"^ is frequency indepen- 
dent, and hence to this order, the single particle spectral function is a simple delta function paik, v\ = 5{v + — Ej^ ). 
Thus G'-^-' has a sharp pole with no lifetime effects. The equations are non linear since the x's are given by 

analogous to Ea (|118p . but with energies defined below in Ea (|120p . The Hartree energies contain a temperature 
dependent renormalization of the band width, via the spin and charge correlation functions. These nearest neighbour 
charge and spin correlations determine ^o, and lead to a shrinking of the band width in case of antiferromagnetic 
correlations. For ferromagnetic, i.e. Nagaoka Thouless type correlations, one has the opposite effect and magnetism 
promotes kinetic motion. Near half filling the density dependent term in Ea (jl20l) is suppressed while the spin term 
survives. Due to antiferromagnetic correlations expected at short distances, the parameter is potentially smaller 
than unity, temperature dependent and can vanish, giving rise to a metal insulator transition The self consistent 
solution of xniq] determines the spin and charge correlation functions. 



D. Hartree Fock Approximation (HFA) 



We next outline the Hartree Fock approximation where the Greens function contains the exchange term obtained 
from Eq (|51|54p by setting the vertices F^ ^> 1 and Fj — > 1. The bare vertices are frequency independent, and all 
vertex corrections are functions of frequency that vanish at high frequencies. Hence the Hartree Fock energy is the 
best possible frequency independent approximation to the correct excitation energy. For this reason, it also gives the 
correct form of the first moment of the Greens function. We write 

G-^Uk] = + m) - Ei"^^ , ^) = i?f ) - i Chf Jk, Chf E cos(g.)m, = {Xfx",n, 

(123) 

where we dropped a k independent term in the self energy (from Sqinq), and assumed the nearest neighbour hopping 
on a simple cubic lattice to simplify the expressions. The vertex functions can be written down from inspection as 

9 

Ft[Pl,P2] = l + ]^^Jpi~qXt[q,q+P2~Pi\ 
q 

rr[pi,P2] = {'^-'^){£pi-£p2) + ^{^q+pi-Sq+P2)^[q]-\'^Jpi-qxT[q,q + P2-pi\- (124) 

q q 

It is straightforward to verify that the Ward identities Ea ipTj) are satisfied exactly, and so this is a consistent scheme 
as well. Unlike the earlier cases discussed, this approximation generates frequency dependent vertices. The vertices 
Ts, t\pijP2] are now functions of the momenta as well as the frequency difference uip^ —^pi- The susceptibility Xhf[Q] 
can be obtained after solving the vertex functions, the singlet and triplet susceptibilities now differ from each other. 
The triplet susceptibility is enhanced at a finite value of Q, whereas the singlet susceptibility is suppressed at finite 



22 



E. Non Linear Hartree (Fock) Approximation (NLH(F)A) 



The next approximation is obtained from Ea (j79|) . where the terms of O(A^) are isolated. Fourier transforming 
Eg ((75)) we find 

^s[Pl,P2] = 1 - Xs[P2 -pi] (,s[Pl,P2], 

3 

Upi,P2] = Y^Xi=c{24. + 42-Pi}' 

^t\pi,P2] = 1 + Xt[P2 -Pl] ^t[Pl,P2], 



^t[pi,P2] = X.oc {2e;^ + J;,_p, + 24 - 3J*^_p, - 28^} 



with xioc given in Ea ljFlip . The "7" vertex" foUow from the stated rules, and we write 

^T[Pi,P2] = (1 - ^)(epi -epj + ^(eg+pi -eq+p^)iy[q] - A^ xT[P2-pi] £.s[pi,P2]- 



(125) 



(126) 



The partner Greens functions of the vertices Ea (|125p . are taken from the Hartree approximation Ea (|119p . Using an 
important corollary of Eq (|9T1) 



i{ujp^ - UJp2)Xs[P2 - Pl] = xT[P2 - Pl 



(127) 



it is readily seen that the Ward identities Eq (inT|) are satisfied for these Greens functions and vertices, and thereby 
the NLHA is also a conserving scheme. We may add the Fock terms to the self energy (the C term in Eq ()123|) as well 
as vertices (the J dependent terms on the RHS of the Ea (|124p ). exactly as we did in going from the Hartree to the 
Hartree Fock theory. This produces the Fock generalization of the NLHA, i.e the NLHFA. We will content ourselves 
with a few comments about the structure of the NLHA below. 

As in the case of the Hartree Fock approximation, the vertices are now frequency dependent. We begin by multi- 
plying Eqs (|125p with G'[pi]G'[p2] and integrating on one of the two momenta to get the susceptibilities: 



Xs,t[Q] ^^GM TsAq-.q + Q] G[q + Q]. 



(128) 



The answer is given as: 



ixs[Q])-' = ixo[Q])-' + j^\'xioc{Jh + 2^%} 

iXtlQ])-' = (xo[Q])-'-^A2x,oc{2J,5-34 + J^ + 2(E^-E^)}. 



using the energy type variables 



E? 



Q] 



Xol 



(129) 



(130) 



and similarly for the triplet channel. A rough approximation is to ignore the frequency dependence of the energy type 
variable E*, E* and to think of them as further renormalized versions of the band energies and the bare exchange. We 
denote these energies with bold letters, to emphasize this point, and to distinguish these from the previously defined 
energies that are all real. Let us first recall from their definitions Eq (|116ll'Tip that both the bubble xo and the local 
susceptibility Xioc are negative variables. If we treat Eg oc eg,, we see that the content of Eqs (|129p is to enhance 
the physical susceptibility (— Xt) and to decrease the physical charge compressibility (— Xs)j by amounts that are 
sensitively dependent on the prefactors. The theory has some resemblance to the random phase approximation, but 
with several coefficients including Xg , Xt that are found self consistently. The results will be published separatelj^. 



F. Frequency dependent self energy 



From the above discussion, in all the schemes discussed so far the self energy involved in the calculations is frequency 
independent. This frequency dependence is important since it provides a measure of the decay of quasiparticles. We 
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see however that the vertices are frequency independent only for the first few approximations, and become frequency 
dependent in the HF, NLHA and NLHFA schemes. Since the self energy ultimately derives its frequency dependence 
from the vertex as in Eq(l5Tl [54l) , we must find a scheme where both vertices and self energy are frequency dependent. 
This is an important problem to be addressed in future work. However, as a via media solution, we may take 
the frequency dependent vertices and substitute them into Eg (1541) . as a non self consistent first approximation. To 
lowest order we find terms that are reminiscent of those encountered in Fermi liquids, with S oc GGG, i.e. with the 
convolution of a bubble susceptibility with a Greens function. This integral is a familiar one from Fermi liquid theory, 
and the imaginary part of the self energy is generically oc oj^. The FS volume is unchanged from the Hartree value, 
in a manner that is quite similar to the standard FL theory. There are several detailed question that remain to be 
worked out regarding the shrinking of the band width and of its temperature dependence, we will return to these in 
futurcSi 



VIII. ANISOTROPIC D-WAVE PAIRING AND SUPERCONDUCTING INSTABILITY OF THE ECQL 

We next study the possibility of a superconducting instability of the ECQL. We study the t-J model without any 
added terms such as phonons, and hence we are looking at the possibility of a spontaneous instability of the type that 
'^He undergoes, when it becomes a superfluid^"— . 

Near half filling, a full Gorkov- Nambu type calculation with anomalous Greens functions within this formalism 
is feasible for this purpose at the level of the NLHA, and we will present the details later^^. To find the existence 
of a d-wave instability, we can take a shortcut; following the precedent in '^He, we extract the effective interaction 
from the NLHA vertices Eq (|125p . Onsite s-wave pairing is excluded by the prohibition of double occupancy, and is 
automatic in the full Gorkov Nambu type scheme for the ECQL, within our formalism^S. However, in the present 
phenomenological scheme, we must implement it by ignoring all but the d-wave channel^. 

The singlet and triplet vertices in Ea (|125l) are at the (NL) Hartree level, and from these we can extract the 
irreducible interactions^^ ^I{p,p';Q) and '^I{p,p';Q) in the two channels, and using crossing symmetry deduce the 
pairing interaction in the particle particle channel. As Leggett points out^, this is achieved more transparently 
when we fit these interactions to a pseudopotential Upp treated at the same (Hartree) level^, and then consider the 
pairing of the pseudopotential Hamiltonian H = Hq + H'^^. A strong short ranged repulsion is added to eliminate 
s-wave channel, this is necessary since we consider singlet pairing unlike ^iJe, where the triplet pairing forbids the 
s-wave channel, by symmetry. The total pseudopotential then is J7pp(l,2) — (oo) x (5(1,2) + [/pp(l,2), and the (spin 
dependent) rotationally invariant pseudopotential given by^ 

C/pp(l,2) = F(l,2)+ai -^2 W{1,2) 

Pl+P'2=P3+Pi 

{PlCri,P2<J2\Upp\p3a3,P4a4) = SaiaJa2a4{{Pl,P2\V\p3,P4) +(71(72 {Pl , P2\W\p3, P4)} 

1(T3 (5(72 <T4(5(Ti(72 {pi,P2\W\p3,P4), (131) 

where c are regarded as the Hartree quasiparticles, and we proceed to find the vertex correction for this Hamiltonian 
at the Hartree level 

rsbl,P2] = 1 + 2 iP2P'l\V\Plp'2) Xs[p'l,p'2] 

Pl+P'2=P2+P'i 

rt[pi,P2] = l + 2 (P2P[\W\pip'2) Xt[p'l,P2l (132) 

Pl+P'2=P2+p'i 

where Xs[p'i,P2] — G\p'i^ s\p'\-,p'-2\G\i[l-^ etc. Comparing with Eq (jl25p we see that 

2(p2,p'l|l^bl,p'2) = -\^is{piM 

2{p2,p[\W\pi,p'2) = +X^^t[pi,P2]. (133) 

We can now insert these potentials into the standard anisotropic Cooper pairing problerr^i^. For the case of singlet 
pairing, we require the pseudopotential for the process \p2 t ~P2 -I) \pi t ^Pi i)- We write the required matrix 
element oi V — W 

1 Tl 

Uppipi,P2) = 22 (1 _„)(1 - s)2 {2 "^0 ^* + (^'■- + '^^t)i^pi + £P2) + Jpi-p^} ■ (134) 
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In obtaining this expression, we have used the definitions of and ^4 in Ea (jl25p . and symmetrized these in pi^p-^. 
It is easy to see that only the third term survives in the d-wave channel, and thus we truncate further 



u: 



pp 



1 



Xs J D 



32 (l-n)(l-f)2 "Pi-P^- 



(135) 



Clearly the same expression holds for the full pseudopotential U' since the strong repulsion has s-wave type symmetry. 

The term XtJp^-p^ cancels out in taking +Ct above, and we observe that the exchange energy Jpi_p2 is multiplied 
by the nearest neighbour density density correlations function Xg = {{npinpi+ij) — "n?) / (n — n^) defined in eg ((56)) . This 
object is closely connected to the pair distribution function discussed in the electron gas, and is well known to have a 
correlation hole, i.e. particles avoid getting close to each other regardless of their spin. In the t-J model, a similar 
depletion is expected so that we expect Xs < 0. If we consider a fully spin polarized liquid, then we can compute Xg 
easily from a Fermi gas picture, and we find Xs = — l(coCjj)|^- This object is negative and small near half filling 
oc —(1 — n). We expect Xs to be negative in the ECQL in the paramagnetic limit, although the magnitude should 
be larger than that for the ferromagnet, since particles need to be neighbours in order to benefit from the exchange 
interaction. 

We now treat U^p^'^"" within the pairing scheme^ and write down the gap equation for d-wave singlet supercon- 
ductivity 



P2 



pip) = 



tanh(i/3£:(p)) 



(136) 



where E{p) = [^^ + |A(p)|2]i/2^ with the Hartree energies — ^0 (cfe ~ A*o)j with /io the chemical potential from the 
first approximation (without the ^0 correction), and 



Co 



n -> ^ 1 



(137) 



as defined in the NL Hartree theory Eq (|120p . This equation can be linearized near the transition temperature, by 
setting E(jp) — ^ |^p| in the summand of the above equation. We assume A(p) oc cos{px) — cos(py), and a simple 
analysis gives the condition for the transition temperature Tc as 



^(cosp^r - cos pyf 



w 



dep{e)ili{e) 



tanh 




tanh 


!Co(ep Mo) 




!Co(ep - Mo) 



1 

a J' 
1 

a J' 



(138) 



where 



V'(£) = -^^{cospx- co&pyf5{ep- e), and 

i-Xs). 



32(l-n)(l-f)2 



(139) 



The density of states p{e) and the angular average ■(/'(e) are easily found for the square lattice in terms of the convenient 
variable u — with W — 4t, and the elliptic integrals E{m), K{m) (where m is the parameter of the elliptic integrals) 



P{e) 



2ttH 



K{1 - 



(140) 



At low temperatures, the sum diverges logarithmically from the region e ^ po- We can extract the divergence by 
expanding the integrand around po, which may be safely taken to its zero temperature limit. We thus find 



1 

To' 



pIpoMpo) log i.upMw^ - pI) 



1 

a J' 



(141) 
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and hence 

keT, ^ 1.13^0 ^J{W^-f^l) e-^o/^* , (142) 

where 

C=a Jpi^^oM^^o). (143) 

This expression is vaUd provided the resulting Tc is much smaUer than the band width 2W, and further we need the 
positivity of the two variables and —Xs- The maximum Tc this approach can yield is ksTmax ~ 1.13^* [W^ — /Iq). 
Taking standard values for parameters in High Tc systems, namely t — 6000*^ and J = 1500°, this maximum Tc 
decreases from lOS^K at n = .75 to 43°K at n = .9, provided we use the Hartree estimates for Xg and with xt = —AAn 
to fit the known ground state energy of the Heisenberg antiferromagnet. However the Tc from Eq (|142|) is much smaller 
than these values, because ^o/^* ^ 1 in the entire range. 

It is interesting to compare our pairing equation (|136I138I142[) with corresponding equations in the work of Baskaran, 
Zou and Anderson (BZA)'^'', who first proposed that the superexchange interaction could lead to superconductivity 
in the t-J model, and to the work of Kotliar— who generalized BZA to the case of d-wave symmetry. BZA's pairing 
equation is obtained from an intuitive argument where the exchange energy is written in a particular factorized way. 
Its mean field theory results in a pairing Hamiltonian that has a striking resemblance to our pairing term U^^^'^'^ . In 
fact their mean field theory transforms to precisely to the above equations if we make the following mappings from 
our calculation: ~^ (1 ~ n), a — 1 and finally adjust the chemical potential /xq — > /xi, with /^i chosen so that 

On the one hand the qualitative conclusions of the two approaches are very close. Within our theory, at least 
within the NL Hartree approximation, superconductivity is possible in the d-wave channel thanks to the sign of the 
correlator Xs\ it turns exchange into an attractive interaction from its initially repulsive character. The mean field 
theory of BZA obtains the attraction by a specific factorization of the exchange energy, and while it is not clear that 
this factorization is unique, it is consistent with the NLHA. 

On the other hand, BZA attain a much greater Tc than our calculation does. This can be tracked down to one slightly 
unfavourable and one crippling difference. The ratio ^/(l — n) is 0(1) and does not make any difference, however the 
ratio p{fj,i)tp{fj,i) / p{^q)^P{^q) is 0(3) for most of the range of densities, and this enhances their Tc somewhat. The 
ratio 1/a is very large 0(30) for most densities, and this makes our Tc come out very small. We thus see that 
the dimensionless constants in our expression Ea (|139p (e.g the factor of 32) make all the difference between the two 
approaches. Our approach systematically leads to these constants as stated, at least within the NLHA, and cannot 
be ignored. We thus feel that superconductivity within the ECQL is very subtle, its currently precarious scale could 
well be influenced by correction terms beyond the NLHA considered here, and must await further investigations. 



IX. PHYSICAL INTERPRETATION OF THE QUASIPARTICLES 



We next discuss the physical significance of the quasiparticles, defined as the poles of G, contrasted with the bare 
particles obtained from Q. The number density of particles n = N/N^ is found by taking the translationally invariant 
limit of 



l + G[i-,i] 

Inverting it, we express the local quasiparticle density ngp 



^{tr G[r,i] -2detG[r,i]} (144) 
^^t*"''^ when V^O. (145) 



nQp[i] = tr Gaa[i ,«], 
n 

(146) 



1 - - 

^ 2 



where we have taken the paramagnetic and uniform limit in the last line. Thus the QP density is always larger than 
the bare density by a factor that is unity at very low fillings and approaches 2 near half (bare) filling. In the case of 
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a general spin population with the partial densities denoted by = Na-/Ns, it is easy to see that 



1 



(147) 



illustrating the fact that both spin populations of the bare electron contribute to that of the quasiparticles of a given 
spin. We see that the the quasiparticle densities come closer to the bare ones as we polarize the t-J model. This is 
natural since the Pauli principle already keeps like spins apart so that the effect of the projection operators is reduced. 
In the fully polarized sector the problem reduces to that of a spinless ideal Fermi gas. 

Fourier transforms can be performed on turning off the sources, since translation invariance is restored, and we 
can construct the occupation in momentum space (see Eq (|113p ) as for the standard Fermi liquid. This is carried 
out in detail in Eq (|114p . In the present case, it is clear from Eq (lll4[ Il47p that the number of electrons "contained" 
in G are greater than those in f/ by a factor -. ^„ . This is the "renormalized" particle number sum rule mentioned 

in Eq (lll4p . From this relation we expect that the other criteria for determining the Fermi surface outlined above 
(after Eq (lll4p ) are similarly scaled. Detailed calculations within various possible schemes are underway, and we will 
comment here on the basis of simple calculations. Within the Hartrce or the Hartree Fock approximations, the self 
energy is real and the various criteria give the same result. The analog of the Luttinger Ward's FS volume theorem3S 
for the extremely correlated quantum liquid (ECQL) holds, provided we replace the electron density by an enhanced 
value as in Eq (|147p . Thus we predict Fermi surface volume H,^^ for the ECQL state, in comparison to the FL (Fermi 
liquid) state and the ferromagnetic (FM) state (or equivalently spinless particles) to be given by 

pg _ n _ n pg ^ECQL _ 2 

"ECQL — 2~n' ~ 2' ~ j-j_Fs — 2-n' 

FL 

This renormalization of the volume by ^* signifies a lack of adiabatic continuity with the non interacting electron 
problem^Sj a key feature of the FL. In Appendix |D1 we locate the origin of the breakdown of continuity. At least 
within the limited setting of the atomic limit of the Hubbard model, we can trace the origin of this change in 
volume. We study the change in functional dependence of the Greens function and self energy upon cranking up 
the interaction strength [/ at a fixed frequency Wn- The distinction between two high frequency limits: the weakly 
correlated (w ^ oo and ^ — 0) or the extremely correlated (oj — ;> cxd and ^ — oo) is responsible for the changed 
volume of the Fermi surface. With the insight gleaned from this exercise, we conjecture the behaviour of the general 
Hubbard model self energy in the limit of extreme correlations (EC). Assuming this conjectured behaviour, we provide 
a variation of the standard arguments^SiSi that yields the renormalized quasiparticle FS volume as in Eq ()148|) for the 
Hubbard model in the t-J regime of parameters. These volumes satisfy the bound 

> ^^CQL > ^fI, (149) 

so that the ECQL Fermi volume differs from both the standard cases for general filling, and approaches that of the 
FL and the FM states at low (n ^ 0) and high densities (n 1) respectively. 

Independently of our proposal, unbiased numerical methods have recently suggested that the Fermi volume of the 
2-d t-J model differs from that of the FL by different enlargement factors ^* as in Eq ()148p . although the factors 
seem a bit smaller—. One curious consequence follows for the nearest neighbour hopping bipartite lattices, e.g. the 
2-d square lattice or the 3-d square lattice. At precisely n = ^, the quasiparticle density is exactly one half. Thus the 
QP FS volume is half of the first Brillouin Zone, and hence they occupy the nested diamond shaped region expected 
for bare electrons at half filling. Beyond this filling, the curvature of the FS changes from electron like to hole like. 
Therefore one would expect the Hall constant of the t-J model to change sign and become hole like at n = |. Studies 
of the t-J model Hall constant^ are consistent with this expectation, showing a change of sign at exactly this filling. 

It is instructive to deduce from the quasiparticle Greens function the time dependent number density ngpfi] = 
tr Gctct[j~,«] and the spin density SQp[i] = ^tr r G[i~ in terms of the bare number density n[i] Eq ()144|) and bare 
spin density S[i] = ^tr f i], with r the three Pauli matrices. These follow from the inverse relation 

^[-''-] = det(l-V,^]) (^-g[-"'-])-g[-"'-]' 



and hence 



nQp\i] = 2 



(l-nH)(2-nH) 
(1 - inW)2 - . 
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Thus the quasi-particle number density is locally related to the bare particle number density and the bare spin density 
in a non linear fashion. These relations are easily inverted as well, 



n t\ = 



2 {2 + (1 - nQp[i]) 7iQp[z]} - J(2 - nQp\i]) [2 - nQp[z] - (2 + nQp\i]) \SQp\i]\ 



5W = Sqp[i] ^ (152) 

The inverse relations are in some sense more fundamental, since the quasiparticles are the basic objects that drive 
the bare particle response. We observe that at a given number density n[z], a tendency to form a local moment by 
the bare particles, i.e. |>S'[z]| ^1/2 n[i\, enhances (reduces) the quasiparticle spin (number) density. Further we see 
that the quasiparticle spin density is scaled down from the true spin density by a factor ^ S = I — n corresponding to 
the hole density measured from half filling, and also the quasiparticle band becomes full when the bare particle band 
becomes half filled. 

It follows from Eq (|151[) that the delicately structured relation ship between quasiparticle charge and spin density 
should be best seen when we dope the uniform ECQL with a charge or spin impurity. In this case the Friedel trapping 
of a single bare charge ends up capturing a non trivial (density dependent) number of quasiparticles, and the spin 
density reflects the charge density as well. This leads us to expect that the role of the impurities would be important 
in revealing the nature of the quasiparticles. One interesting feature is that when the bare particle density is close 
to unity at any point, the quasiparticle spin density has a local minimum at that point, and thus displays a non 
monotonic behaviour. 

From the above construction, we conclude that the charge of the quasiparticles qqp must be regarded as a density 
dependent fraction of the bare charge qe 

qqp = {1 - ^n} (153) 

in order that the total charge remain invariant; i.e. qqp Nqp — qe N. Near half filling, the charge of the quasiparticles 
is ^ i. Therefore the flow of a bare particle is equivalent to that of a sufficient number (23;^) of quasiparticles, so 
that the total charge is balanced. These fractionally charged particles are defined in the many body context without 
any specific single particle basis. These fundamentally arise in terms of a modified Pauli principle implied in the 
equations Eqs ()114ll46|147|) . 



X. CONCLUSIONS AND SUMMARY 



In this work, we have presented a systematic study of the t-J model by using the Schwinger technique of source 
fields. In addition we have developed a specific methodology to overcome the problem of non canonical Fermions 
forced upon us by the infinite U constraint in the model'^^'S. Since the method is technically quite involved, we have 
presented the details in a self contained fashion. We obtained the exact Schwinger Dyson equation for the t-J model, 
and hence a closed expression for the inverse Greens function G^^[k] in terms of the vertex functions. Both singlet 
and triplet particle hole vertices are needed to complete the definition of the G~^[k]. The vertices are reported up to 
the neglect of the source derivatives of the vertex- this is a natural stopping point since we need to first understand 
the consequences of the many terms generated so far. In order to facilitate concrete approximations, we presented 
the Ward identities for the current and density vertices. 

The resulting Greens functions and vertices form a hierarchy; this is in many respects similar to the one usually 
encountered for standard Fermi systems, but with extra features that arise from the dynamical analogs of Gutzwiller's 
projection operators^ for the ground state. As with standard models of weakly interacting Fermions, such a hierarchy 
is the proper setting for exploring other features that might lead to controlled approximations. The existence of a 
low density nuclear matter (Brueckner, Galitskii, Migdal) type or high density electron gas (Bohm, Pines, Brueckner, 
Gell-Mann) type approximations in the Fermi systems are examples of such an emergent process, and are described 
in various texts^'^'^'^. Our preliminary search shows a natural ordering of terms in the vertex, where the hole 
density plays a central role. This scheme is currently under numerical study and results will be reported later---. 

The statistical mechanical equilibrium state underlying our Greens functions is the extremely correlated electron 
liquid. This quantum liquid breaks no spatial or temporal symmetries. It has specific signatures that distinguish it 
from the Fermi liquid. In particular, the Fermi surface volume naturally differs from that of the Fermi liquid. The 
elementary excitations of this liquid are best viewed as fractionally charged quasiparticles, whose charge is determined 
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by the density. At all densities of the t-J model, the particle number sum rule requires a Fermi surface that is larger 
than the standard Luttinger-Ward Fermi surface22., by an amount ^* = t^-tt- In Appendix |Dl we have presented a 

2 

suggestive variation of the standard argument for the Fermi surface volume that gives us the scaled volume. Further 
studies are needed to compute the spectral functions and thereby answer the issue of sharpness of the quasiparticles in 
a complete fashion, by computing G~^[k] within a controlled scheme. The fractionally charged quasiparticle picture 
of the ECQL is a quantitative description of the so called Lower Hubbard BanS^.. 

An interesting feature of the theory is that the (^"^[fc] involves the static density and spin correlations at nearest 
neighbour separation; this leads to the narrowing of the bands and the possibility of a metal insulator transition 
near half filling driven by local antiferromagnetic correlations (rather than true antiferromagnetic order) . In general 
the ECQL is prone to various instabilities such as the antiferromagnetic, insulating or superconducting states. We 
presented a calculation of the superconducting instability towards d — wave order in Section I Villi the final formulae 
have considerable similarity to the RVB theory^!, but end up with a much lower due to some dimensionless factors 
that are unavoidable in this theory. More generally, in the ECQL, unlike the standard Fermi liquid instabilities, one 
does not need to deal with a large energy scale such as U, since the Hubbard operators already deal with the local 
constraint efficiently. This is a great advantage since in the FL, the energy scale U skews the picture of instabilities 
by overemphasizing the magnetic instabilities. 
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Appendix A: Identities involving the source derivatives. 

Let us next work out the transformation of the derivative 



A-^H • D[j] • A[fc] - A-'[i] ■ D[j] ■ A[k] + -.A-'li] ■ D[j] ■ A[k]: (Al) 

where the first term consists of D[j], the overline symbol is analogous to a contraction in field theory. In the present 
context it implies that D is acting as both a matrix as well as a derivative on A[fc]. In the second term the vertical 
dots denotes "normal ordering" w.r.t. the derivative operator, i.e. in matrix element form, we take the derivative 
operators to the right of the A[k]. These two terms arise from the action of Z? as a matrix derivative acting upon the 
A[k] term, and also on whatever stands to the right of the expression, and is the analog of the covariant derivative in 
non abelian Gauge theory. We will next show that this may be expressed as a useful identity: 

A-i W • m ■ A[k] - ^l[^,J] ■ k] + ■.^l[^,J] ■ V[j] ■ ^[j, k]\, (A2) 

where is defined below in Eq (|A3p . it is the transformed version of D in terms of the transformed source terml^. 
Let us first prove an identity for the derivative. 

■.A-\{\ ■ D[i\ ■ A[i\: = V[i\ 

This useful identity implies that in calculations where the source l/l is turned off, as for e.g. in Appendix IF] where we 
calculate the susceptibility relations, we can ignore the distinction between D and V since A 1. 

We use Eq ((24l) i.e. A^^[i] = jc/^^j] A'^li] to rewrite the inverse A. Taking components, we write this equation as 

A-^W ■ ■ A[z];) = a„a,A^^,J.]A.,.,[z] (A4) 

Since Vi — A^^[i] • Vi • A[i], we may write 



S 1 fiA n ^ 

5V^ " detAW W^-.-.W^ 



^ -A.„.Jz]A^,„Jz]— 1^, (A5) 



detAfil '^"'^'^^^ 'Td<^, L J ^^ff^ffd : 
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where we have using the definition of the conjugate in going to the second hne. Substituting into Ea (jA4[) we see that 
the spin components are now properly arranged to yield delta functions so that the identity Eq (|A3p is proved. 
From the identity Eq ()A3| , we may write another useful transformation: 

■.D[i\':^':A[i\-V[i\-A-^[i]:, (A6) 

this helps us to replace functional derivatives w.r.t. the original source V with those w.r.t. the transformed source V. 
Using this, we may at once rewrite the LHS of Eq (jA2[) 



;A-iW • D[j] ■ A[ky. = -.A-^ii] ■ A[j] ■ V[j] ■ A-^[j] ■ A[k]-. 

= ■.ti[t,j]-V[j]-ti[j,k][ (A7) 

This proves the second term in Ea (|A2[) . 



1. The I' matrix. 

We turn to the first term in Eq (jA2p . we see from Eq([6T|) that 

7F] 



.[fc] = -^(l-G'=[A;-,fc]), (A8) 
where we denote 7[fc] = 1 — det G[k^ , k], therefore 



iy[k,j] = A-'[j] ■ D[j] ■ A[k] = A-'[j] ■ :A[j] • V[j] ■ A-^j]: ■ A[k] 

= ■.V[j]-A-^[j]:-A[k] (A9) 

In component form we write, 



''<T^,a2[kJ] = A„^,,^[j]V^,^^^[j]A^,^^^[k], 



l[k] 



(AlO) 



where we used 



^ detG[A:-,fc]=Xar.'Jfc,fc;j] (All) 



For spin diagonal sources, we find 

y..[k,]] = ^ [A'M ((Aaa[fc]G..[fc,fc] - l) x^'\k,k,]]+A„,[k]G,^[k,k]x^^\k,k,]]) - A^l[j]x^''[fc,fc,j1 



l[k] 

Upon turning off the sources, we find 

1 - n ''2 ^'^ ' 2 



(A12) 



i^a,^.M,b] = <5,,,,, ^{(5-l)x(i)[a,a;6] + 5 X^'^Ka;6]-x('Ha,a;&]}- (A13) 



2. The e matrix. 

We calculate the object 



e[r,.s,m] = \V[r]-fi[r,s]\ ■ G[sJ].G-'[f,m], (A14) 
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this appears above. Taking the matrix element of this we write 

= aiaalJ.a^,c7dr,s]Gai^aAs,p]^9l%l\p,m;r]. (A15) 

For spin diagonal sources we find 

eaAr,s,m] = A,,[s]A;,i[r]Ga,[s,A:]r(2)[fc,m;r] - As5[s]Ar-i[r]Gs5[s,fc]r(3)[A:,m;r]. (Af6) 

Upon turning off the sources, this becomes 

ea,,aAr,s,m] ^ S,,,^,G[s,k]{r^^\k,m;r]-r^^^k,m;r]} 

= 5,,^,,G[s,A:]r(P)[fc,m;r]. (Af7) 

Appendix B: The $i matrix. 

We study the properties of the matrix <i>i defined in Eg ([5^ as 

= A-i[z].a,.AH, (Bf) 

and begin by starting from Ea ipO]) 

a,,r(e-^Xf-"^(T,)) = -T(e--^[Xf"'"'(r,),i/]) +Vr"^(T.)r(e-^[Xf^^=(r,),Xf"''^^(r,)]) and hence 

{{{X^r\n)X^"^{n))) ((Xf"0(r,)X^(r,))) ) 

+\ 'hi mr'in)xr'in))) {{Xl-''^{n)xr^{n))) ) (B2) 

Thus we obtain on using the definition /S.ai,a2[i] ~ ^(Ti.a2 ~ {(^1(^2) {{^1^'^^ {'''i))) ■, a-nd the k conjugate of the source 

{a,a,) {{{xp"in)xf^^^{n))) {{xr'{n)xr^{n))) ) 
-\ J^,i (^1^2) {{{xr^-{n)x!^^'Hn))) - {{xr^{n)Xf^-^Hn))) ) , (B3) 

and after some calculation we find 

= Ff^-VfVi. (B4) 
= t[i, /] A-i W • (^^'^ I] - G'-"^ [I, *]) + ^Ah l]{m - I] ■ AW) 

+ y[i,l]A-'[^-{D[^-A[l]-D[l]-A[2]). (B5) 

The second and third lines of Eq (IB3p correspond to the term $ in Ea (jB5[) . and it should be noted that the time 
arguments of the X's are all equal. From this it is clear that the term $ and vanishes on turning off the sources, due 
to the cancellation of terms obtained by interchanging i -o- 1^. 

Appendix C: Rotation Invariance and the Nozieres relations. 

We summarize the rotational invariance argument for the various vertex functions, these are in close parallel to 
the standard Fermi liquid arguments as related by Nozieres^. We further assume that time reversal and parity are 
also preserved in this putative ECQL state. Let us consider the case of an extremely correlated quantum liquid with 
no broken symmetries. Since there is no source term, rotation invariance leads to a useful spin decomposition of 
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the vertices. A general four legged object XaaaJ" analyzed initially, and we show that other objects have the same 
decomposition. Since the symmetries of this hquid are identical to those of the Fermi liquid, we can extract the 
decomposition by studying the parallel problem. We schematically view this object as the scattering amplitude for a 
particle and a hole, from the definition of this object, it is related to the correlator (assuming a suitable set of particle 
times) 



(CI) 



where /j^ creates a particle with spin cr, and h\ = a creates a hole with a spin a. The hole creation operator has 
by standard convention an extra factor of cr, in order to ensure that the spin flip operator of the particle f^J^ maps to 
the spin flip of the hole /ij/i^, rather than its negative. This object is Fourier transformed according to the standard 
rules Eq (|Hip . and we will omit the momentum labels for brevity. 

We may thus view the x as a scattering amplitude for a process taking a particle hole pair with a certain initial 
state to a final state^'^ as follows: 



Initial state = {cd, Cc} 
Final state — {<yb,Oa}. 



(C2) 



We first note that due to time reversal invariance, the object satisfies 



(C3) 



with reversed momenta. However, using parity, we can reverse all momenta again restoring them to their original 
values, therefore this relation is true with fixed momenta. 

Let us call S^*°* as the total spin of the particle hole pair. We must conserve the z component of this object 

We must also conserve the total spin magnitude \S*°*\ — or 1 in this process, since the ground state is rotationally 
a singlet. For implementing this, we further decompose the scattering into the particle hole singlet and triplet channels 
as follows. The particle hole states can be represented in one of two possible schemes that are illustrated for two of 
the states as follows: 



Scheme A 


Scheme B 


Singlet m-it} {f\hl-flh\} 


{f\h + flfi} [tt + U] 


Triplet [U + m {flhl + flh\} 


iflh-flh) [tt-U] 



(C4) 



Using Scheme A, we write the four possible states (displaying the initial state) and the corresponding scattering 
amplitudes as 





Singlet Channel 

2 X/CTcr' Xa'a' ~^ 




{tn or {;;} 


Triplet Channel 

Xt = x%% ^ 





(C5) 



where the non zero amplitudes are denoted as 



^'^-Uxs + xt) 



[Xs - Xt) 



X^'^=Xf 



(C6) 
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Assuming spin isotropy of the quantum liquid phase, and all states of a given spin must give the same results for the 
scattering amplitude. From the two possible states of the triplet channel, we glean an important and labour saving 
identity: 

We find in the equations above that a specific combination for the susceptibility x the vertex T arises repeatedly: 

A AfTcr AfTcr A A •> 

= \{Xs-ixt}. (C8) 

This corresponds to the correlator of singlet Cooper pairs oj particles, i.e. q]6[r, s]), with g] ~ f^{p)fl{Q) — 

so we call it the pair channel. 
We also see that the vertex function F "inherits" this symmetry Ea (jC7p of x- This follows from the definition in 
Ea PO]) and the isotropy and diagonal nature in spin indices of the G's i.e. Ga-ia^ 5^^,(7 jG^^ on switching off the 
source terms. 

Appendix D: The Hubbard model atomic limit and FS volume 

In this appendix, we locate the point where adiabatic continuity is lost, in going from the Hubbard model to the 
t-J model. The context of the discussion is the atomic limit, where we can trace this quite explicitly. However, the 
generality of the argument regarding the high frequency limit and the asymptotic form of the Greens function and self 
energy in the two models is to be noted. One of the consequences is that we have identified the essential problem in 
applying the Luttinger Ward theorem to the t-J model, it must suffer corrections due to the neglect of the boundary 
term n2 as outlined below in Ea (jD9|) . 

We recall the solution of the Hubbard model in the atomic limili^ t ^ where 



Gatomic(j<^n) 



1 - - 
2 



- H iuJn +pi—U 
1 



with 



= u'l + U'-^^^-fZ—ry (Dl) 

2 lUJn + ^I,-U(l- ^) 

Here the symbols G, S stand for the standard definitions of the Greens function and self energy^^'^'* and should not be 
confused with the t-J model objects defined in this work'^'*. We are interested in understanding how these functions 
evolve, as we go from weakly correlated (WC) to extremely correlated (EC) limits. In the EC limit we recover the t-J 
model solution, of course in the trivial limit of i, J — > 0. For this we need to understand two distinct high frequency 
(HE) limits 

WCHF hmit LOn-^oo f/ - 0(1) — ^ 

ECHF limit w„^oo f/->oo — -s-oo. (D2) 

We observe that in the WCHF limit, G ^ l/itOn and E ^ cl + c2/ia;„, and thus behave exactly as one expects^i^. 
The canonical nature of the anticommutation relations of the Fermions fixes the coefficient of 1/ itUn in G as unity. In 
the EC limit, we send U ^ oo first, so that 

1 - - 

Gatomic-Ec(«'^n) = ; > (D3) 

lUJn + M 

the coefficient of l/iojn is now 1 — § expected behaviour from Eg (1^51) . this refiects the non canonical nature of the 
electrons in this limit. This coefficient is a statement of the density dependence of the number of states of the "lower 
Hubbard band" . The self energy in the EC limit Ea (jD2[) is easily found to be 

^EC 



{iujn) = co(iw„ + m): where cq = . (D4) 
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Reflection shows that this feature of a hnearly growing E as a function of ia;„, an initiaUy surprising resuh, is natural 
in this case. A hnear growth of E with iw„ with the exact coefficient cq in Eq (|D4p . is needed to match the asymptotic 
high frequency behaviour of G in Eq (jD3P with the correct and related coefficient. Thus we see that the loss of 
continuity in functional form between weakly correlated and extremely correlated Greens function G's and S, arises 
due of the order of limits implied by Eq (|D2p . The EC limit corresponds to taking J7 — > oo before any other limit; this 
process tosses out a fraction of the states of the Hilbert space, and thereby redefines the asymptotics of the remaining 
scales. We conjecture that this behaviour of the Hubbard model self energy S is inevitable in the generic case of finite 
hopping as well. Thus our conjecture is that in the extreme correlations limit of the Hubbard model, the self energy 
must have a form (z = iijJn) 

2)/limu^oc <^oiz + m) + SRogular(fc, z), (D5) 

with Co given in Eq (|D4p where ER,oguiar(fc, z) is a well behaved function (i.e.^ ci + C2/2 ) at high z. For consistency, 
it must in fact be related to the self energy S[k, iuJn] of the t-J model in Eq ([54| with Jij — > through the relation 

SRcgular(fc,it^«) = T— IT ■S[k,iu}n]. (D6) 
^ 2 

Armed with the conjectured Eq ljDSp . and another assumption detailed below, we can modify the original argument 
of Luttinger and Ward^S to obtain the volume of the ECQL Fermi surface. Following the masters, we begin with the 
relation between the number of particles and the Greens function (with 77 = 0^) in a FL; 



= 2^G(fc,icj„)e'""'' 



A; 



? y_ — {G(fc,x-f77)-G(A:,a; + iry)}. (D7) 



-''s ^ J —00 
k 



In the first line, the frequency sum is replaced by a contour integral after multiplying by the Fermi function. Next we 
deform the contour to run parallel to the real line, and taking the T=0 limit gives the second expression. From the 
definition G{k, z) = l/{z + ji — Sk — E(A;, z)), we write 



G(k,z) = -4^1ogG(fc,z) + G(fc,z)-^S(fc,z), 
dz dz 



so that n = ni + n2 with 



-1 = -WXL^ £ {^logG(fc,x-^^)-^logG(fc,. + ^,)} 

A; 



1 2 ^j^g G {k, ^ 00 -ir])G{k, if]) 



2m Ns ^ G{k, -00 + i?7)G(fc, -i-q) 

" k 

where Q{x) is the usual Heaviside function (1 or 0). This evaluation is parallel to the one in the original argument, 
wherein rti is the sole contribution to the volume theorem. The number ni, and hence in a FL the total density n, is 
found by adding up the number per volume of k values where the G(fc,0) at the chemical potential is positive. We 
next consider 

71,2 = J ^ |G(fc, a; - i77)-^E(fc,x - ir;) - G(fc,x + ir?)^E(fc,a; + . (D9) 

k 

This term is usually integrated by parts-^-'^^ and the boundary term E(fc,x)G(fc,x) discarded at infinity, assuming 
that the growth of S(fc,x) is slower than linear in x. From our discussion above, this assumption is incorrect for the 
EC or t-J limit. There the linearly growing T,(k,x) precisely catches up with the inverse linearly decaying G{k,x), 
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giving a non trivial boundary contribution. The corrected answer is most easily found if we use our conjecture Ea (jD5[) 
for decomposing T, into two parts, and we find 

ri2 = CQU + n^, with 



k 



(DIO) 



Now 713 can be integrated by parts safely. We further assume with Luttinger and Ward that SRcguiai (^, z) — "^^"S"'" 

with a suitable functional $Rcguiar- This is the second part of our conjecture alluded to above, and is based upon the 
idea that the Luttinger Ward functional ^lw itself can be decomposed into a regular and singular parts'*'*. With this 
it follows that 713 — 0, as shown in^"^. Thus we find in the EC limit 

n = ni + cq n 

TL 

= (1 — 2''''^i' ^^"^ thus 

^ ^e(G(fc,0)), (Dll) 



1-f iV. 



where we used Eq (jD4p in obtaining the second line from the first. We emphasize that this renormalized version of 
the Luttinger Ward volume theorem, is only applicable in the EC limit as in Ea (jD2[) . and of course is the same as 
the FS volume for the ECQL in Eq([I151). 



Appendix E: The different limits of vanishing source terms 

There are many possible ways to turn off the sources. The generic one used for most calculations is Case.l below, 
and we also list some other possibilities for completeness. We may distinguish three possibilities: 

• Case:l All source terms 

where ^ is infinitesimal. In this case we can readily establish that 

The same holds for V — oc (^^) and also D — T> oc (^^). Hence to linear order in ^ we may ignore the distinction 
between V and V. 

• Case:2 Spin diagonal V where 

In this case again we see that V = V since all off diagonal Q^s vanish identically. This turns out to be not 
interesting since the G's remain off diagonal in spin space here, because the total source term is not just V ( 
see next case). 

• Case:3 The spin polarized case where we choose 

but much smaller than V^"^ — V^^, the latter being set to zero last. In this case, the symmetry of the underlying 
ECQL state gives us Ga;9 oc ^, and Ga,a — Qs.s is arbitrary. Thus we allow for ferromagnetic polarization in this 
case. The total source term from eq (|B4|4ip is Vtot = W + y^'^-' — V'-'^''. We compute this to first order in ^ and 
find 



mi 
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The other two components can be read off from the transformation generated by reversing all spins, i.e.t-^-^i- 
From this equation, we learn that the G's will be diagonal in spin space if we adjust 

(1 - gur,z])v^^ = Qltr,^]iv;^ - k^^), (es) 

since the G's respond to the total source term. 

Appendix F: Various Susceptibilities and their relationships and Sum-rules. 

We need the relationship between the response function T of the Greens function Q to the external potential V and 
the function x defined in Eq (p9|) . 
Let us write basic identity: 

ft - - 

= (1 - 1^)X''.:X [P, 9; r] - aaC7tT2% [p, p; r] G.,., [p, q] (Fl) 

where we have used the relation ^7[p, 9] = A[p].G[p, g] and the relations Ea (jA6l) . Therefore specializing to the non 
vanishing spin configurations we find 

T«[i,j;fc]-(l-^)xW[z,j;fc] = -G[t,j]T^^^r,t;k] 

T(2)[z,j-;fc]-(l-^);^(2)[j,^-;fc] = -G[*,j]TW[z-,*;fc] 



T^^\t,j;k]-{l--)x^^\t,j;k] = +G[*,j]T(3)[i-,j;fc] 

and Fourier transforming, (with T[Q] = J2q ^[l^ Q + Q]) 

TW[pi,p2]+T(2)[p2-pi]G[p2] = (l-f)x^'^bi,P2] 
T(^^[pi,P2]+TW[p2-Pi]Gb2] = (l-|)x^'^bi,P2] 

The susceptibilities at finite wave vectors follow by setting P2 = Pi + Q and summing over pi so that 

2 — n 2 

T(^^[Q] + :7^tW[Q] = (l-;)x(^)[Q] 
2 — n 2 

2 — n 2 

so that 



(F2) 



(F3) 



(F4) 



1 1 — n 

Xs[Q] = -^^—jr^^siQ], xt[Q] = _ n^2 ^*['3] 

n (1 — — )^ 

2 1 — n 

(F5) 

In order to gain intuition for these objects, we note that the "physical" (i.e. positive definite) magnetic and charge 
susceptibilities (i.e. compressibility^) at finite wave vectors are given by 

XsptniQ] = -2/x| Tt[Q] 

XchargelQ] = - ^s[Q], 

(F6) 
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with the normahzation that the corresponding objects for the (non interacting) Fermi gas are Xspin = 2/Xgn[0] and 
Xcharge = ^[0] respectively, with n[0] the density of states per spin per site at the chemical potential. Notice that the 
triplet object Xt turns up with an explicit factor 1 — n, and so it vanishes near half filling, and conversely the physical 
magnetic susceptibility is obtained from it by dividing with this factor. 

From the definitions T[i~,i;r] — ^^p-^/\/-i.o (with suitable spin indices), we can relate the susceptibilities in real 
space to physically interesting correlations 

^^%-.r,3\ = -((Xfxf)) + ((xr)) ((4^)) 
T(3)r,z;j] = -((Xfxf)) + ((xf))((xf)) 

(F7) 

We next turn off the sources, thus at equal times = Tj , we may write the correlation functions in terms of the 
physically meaningful charge and spin correlators: 

1 2 ^ ^ 

If we also set i = j, we to obtain the local and equal time susceptibilities 

T« = -(- - 1) t£) = i-r T^i = -- 

Ts,..^--{l~n) T,,,„. = --. (F9) 

These are useful for the sum rules that we discuss next. 

a. Susceptibility sum rules 

The local x's follow from Eq ljFSP by summing over Q, and these provide us with sum rules for the susceptibilities. 
We find that both Xs and Xt satisfy exactly the same sum rule: 

Q 2) Q 



(FIO) 



with 



_ n (1-n) 

X- - -2 j^^. (Fll) 

The vanishing with 1 — n of both the spin and charge x's is an interesting consequence of our construction. The 
equalized local singlet and triplet sum rules simplify further analysis. The two T's are the physical susceptibilities, 
relevant for neutron scattering, NMR, charge response and other probes. These are naturally distinct from each other; 
for example at the insulating limit n = 1 there is a non vanishing spin response^ but no charge response. 

Appendix G: Detailed Vertices 

In this section we present the vertices, where we have dropped the higher order vertices, i.e. set |^ — 0. We first 
break up the self energy given in Ea (j43p into convenient smaller terms, and then present the the singlet and triplet 
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vertex corrections arising from these six terms. 

= -t[i,j] {A[j] - (1 - |)), S2[i,j] = -t[i,j],y[i,j], 

1 Tl 

SsliJ] = -t[i,k] @[i,k,j], S4i,j] = S[i,j]-J[i,k] {fi[i,k] ■ A[i] - (1 - -)}, 

S5[i,j] = + ^J[i,k] iJ.[i,k] ■ i'[k,i], Se[i,j] = ^J[i,k] iJ.[i,k] ■ @[k,i,j]. 

(Gl) 

We pull out the explicit factors of A = and present the answers as a series in A. The answers are given in real 

space with the external space time indices i,j,m and four internal indices (summed over) a,b,c,k. 
We first write the singlet vertices Ts[i,j;m\r = '^'gy^'^^ with 1 < r < 6: 

Ts[i,j;m]i = ^(n - 2)H[i,j]xs[j,j,m] 

Ts[i,j;m]2 = ^{n - 2) A t[i,j] (^{n - 2)xt[i,i,j]Xs\j,j,m] - 2rt[a,b,j] (G[i,a]xs[b,i,m] + G[b,i]xs[i,a,rn]j^ 

+ 2)* [«, i] (-2rs [a, b, j] (G[i, a]xs [b, i,m] + G[b, i]xs [i,a,m]^ - {n- 2)xs j] {^Xs [i,i,m]-Xs [j, j, r. 

3 

+ g - 2)^n X^t[i,j]xt[i, i, j]Xs[i, i, m] 
^s[i,j;rn]3 = ^t[i,k] {rs[c,j,i] -3Tt[c,j,i]) (^{n - 2)G[k,c]{x8[i,i,m] -Xs[k,k,m]) - 2xs[k,c,m]j 
Ts[i,j;m]4, = —{n - 2f5[i,j]J[i,k]xs[k,k,m\ 

^s[i,j;rn\5 = ^(n - 2) A d[i,j]J[i,k] (rt[a,b,i] (G[k,a]xs[b,k,rn] + G[b,k]xs[k,a,rn\^ - {n - 2)xt[k,k,i\xs[i,i,m]^ 

+ ^in- 2)d[2, j]J[i, k] (2Ts[a, 6, i] {G[k, a]xs[b, k, m] + G[b, k]xs[k, a, m])) 

- (n - 2f5[i,j]J[i,k] {xs[k,k,i] {2xs[i,i,m] - 3xs[k, k,m])) 

- A(„_2)2(2n-1) A^ 5[i,j]J[i,k]xt[k,k,i]xs[k,k,m] 

rs[«,i;m]6 = ^J[i,k]xs[i,c,m]{r4c,j,k] - 3rt[c,j,k]) (G2) 
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We next note the triplet vertices 

^t[i,j;m]i = -i(n-2)2 >^ t[i , j]xt[j , j , m] 
rt[i,j;m]2 = |^(n-2) Af[i,j]| x [ 

2G[i,a\xt[b,i,m] {T,[a,b,j] - {n + l)Tt[a,b,j]) 
+ 2G[b,i]xt[i,a,m] {r,[a,b,j] - (n - 3)Tt[a,b,j]) 

- (n - 2) (xs j]Xt\j, j, m] + 2xt[i, i,j]xt [i, i, m])] 

- ^{n-2fX^t[i,j] {nxs[i,i,j]xt[i,i,m] - Xt[i,i,j]Xt\j,j,m]) 
^t[i,j;m]3 = - 2) A G[k,c]t[i,k] {xt[i,i,m\ - Xt[k,k,m\) {Ts[c,j,i]+Tt[c,j,i]) 

- ^t[i,k]xt[k,c,m] (r^fc, + T^c, 
Tt[i,j;m]4 = ^(n- 2)^ A 6[i,j]J[i,k]xt[k,k,m\ 
rt[i,j;m]5 = S^^in - 2) X 5[i,j]J[i,k]j x [ 

G[k, a\xt[b, k,m\ {{n + l)Tt[a, b,i\ - Ts[a,h,i\) 
+ in- 2)xs[k,k,i\xt[i,i,m] 

- G[b, k]xt[k., a, m] (F^a, b, i]-{n- 2,)Tt[a, b, i])] 

- ^(n - 2)^XH[i,j]J[i, k] {xt[k, k, m] ((1 - 2n)xs[k, fc, i] + (2n + l)xt[k, fc, ^]) - 2xt[k, fc, t]xt[i, m]) 

Tt[i,j;'m]6 = ^J[i,k]xt[i,c,m] {rs[c,j,k]+Tt[c,j,k]) - i(n - 2) A G[i,c]J[i,k]rt[c,j,k] {xt[i,i,m] -Xt[k,k,m]) 

(G3) 



39 



Appendix H: Fourier transform convention. 



Our convention for various Fourier transforms is summarized here. 



pr ~ p.r — LUpT 

- hu 

k - UJ„ k 

G[a,b] = ^e^'=("-'') G[k] 

k 

Pl,P2 

f/3 fP 



Pl,P2 

Q 

X[Pl^P2\ = G[pi]T[pi,p2\G[p2\ 

p 

p 

xiQ) = Y.^^P^P + Q) 
p 

T[Q] = Y^[p,p + Q]. (HI) 
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